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Abstract 

In these lectures I present the key ideas driving the field of relativistic heavy- 
ion physics and develop some of the theoretical tools needed for the description 
t^- ■ an d interpretation of heavy-ion collision experiments. 
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1 PROLOGUE : THE BIG BANG AND THE EARLY UNIVERSE 



Matter as we know it, made up from molecules which consist of atoms which consist of electrons circling 
1—5 | around a nucleus which consists of protons and neutrons which themselves are bound states of quarks 

and gluons, has not existed forever. Our universe originated in a "Big Bang" from a state of almost 

infinite energy density and temperature. During the first few microseconds of its life the energy density 
I , in our universe was so high that hadrons (color singlet bound states of quarks, antiquarks and gluons), 

J> | such as the nucleons inside a nucleus, could not form. Instead, the quarks, antiquarks and gluons were 

deconfined and permeated the entire universe in a thermalized state known as quark-gluon plasma 
^ . (QGP). Only when the energy density of the universe dropped below the critical value e cr ~ 1 GeV/fm 3 

and its temperature decreased below T CT « 170 MeV, colored degrees of freedom became confined into 
2 ' color singlet objects of about 1 fm diameter: the first hadrons formed. 

After the universe hadronized, it took another 200 s or so until its temperature dropped below 
~ 100 keV such that small atomic nuclei could form and survive. This is known as primordial nucleo- 
Qh! synthesis. At this point (i.e. after "The First 3 Minutes") the chemical composition of the early universe 

was fixed ("chemical freeze-out"). All unstable hadrons had decayed and all antiparticles had annihi- 
lated, leaving only a small fraction of excess protons, neutrons and electrons, with all surviving neutrons 
bound inside small atomic nuclei. The chemical composition of the universe began to change again only 
. \ several hundred million years later when the cores of the first stars ignited and nuclear fusion processes 

^ ■ set in. 

After primordial nucleosynthesis the universe was still ionized and therefore completely opaque 
to electromagnetic radiation. About 300000 years after the Big Bang, when the temperature had reached 
about 3000 K, electrons and atomic nuclei were finally able to combine into electrically neutral atoms, 
and the universe became transparant. At this point the electromagnetic radiation decoupled, with a per- 
fectly thermal blackbody spectrum of T w 3000 K ("thermal freeze-out"). Due to the continuing 
expansion of our universe this thermal photon spectrum has now been redshifted to a temperature of 
about 2.7 K and turned into the "cosmic microwave background". The number of photons in this mi- 
crowave background is huge (about 250 photons in every cm 3 of the universe), and they carry the bulk 
of the entropy of the universe. The entropy-to-baryon ratio of the universe is S/A ~ 10 9±1 ; its inverse 
provides a measure for the tiny baryon-antibaryon asymmetry of our universe when it hadronized - a still 
incompletely understood small number. 

The only other surviving feature of the Big Bang is the ongoing Hubble expansion of our universe, 
and the structure of its density fluctuations, amplified over eons by the action of gravity and reflected in 
today's distribution of stars, galaxies, galactic and supergalactic clusters, and dark matter. Using these 3 
or 4 observational pillars (today's expansion rate or "Hubble constant", the microwave background spec- 
trum and its fluctuations, the primordial nuclear abundances and, most recently, also today's spectrum 
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of density fluctuations) together with the equations of motion of general relativity, we have been able to 
reconstruct the cosmological evolution of our universe from its origin in the Big Bang. However, try as 
you want, we will never be able to directly see anything that happened before 300 000 years after the 
Big Bang, due to the opacity of the Early Universe. In particular, the all-permeating quark-gluon plasma 
which filled our universe during the first few microseconds will always remain hidden behind the curtain 
of the cosmic microwave background. This is were relativistic heavy-ion collisions come in: It turns 
out that we can recreate this thermalized QGP matter (or at least some decent approximation to it) by 
colliding large nuclei at high energies. To elaborate on this is the subject of these lectures. 

2 A FEW IMPORTANT RESULTS FROM LATTICE QCD 
2.1 Lattice QCD in 3 minutes 

We know from lattice QCD that the quark-gluon plasma exists (see Ref. [1] for a recent review). Lattice 
QCD is a method for calculating equilibrium properties of strongly interacting systems directly from 
the QCD Lagrangian by numerical evaluation of the corresponding path integrals. One starts from the 
vacuum-to-vacuum transition amplitude in the Feynman path integral formulation 

Z = J VAl{x)V^{x)V^{x)e^ dixC[Aa ^ ] , (1) 

where the phase factor depending on the classical action / d 4 xC[A, xp, tp] is integrated over all classical 
field configurations for the gluon fields, A^{x), and quark and antiquark fields, ip( x ) an d 4>(x)- The path 
integral is dominated by those field configurations (A fJ/ (x),ip(x),ip(x)) which minimize the classical 
action and render the phase factor stationary, i.e. which satisfy the classical Euler-Langrange equations 
of motion. These classical solutions define the classical chromodynamic field theory. Dirac and Feyn- 
man showed that integrating over all field configurations instead of only the solutions of the classical 
equations of motion produces the corresponding quantum field theory. The Pauli principle for fermions 
is implemented by postulating that the classical fermion fields ip(x) and tp(x) are Grassmann variables 
satisfying ip(xi)ip(x2) + ip(x2)ip(xi) = etc. 

Starting from Eq. (1), one obtains an expression for the grand canonical partition function of an 
ensemble of quarks, antiquarks and gluons in thermal equilibrium by an almost trivial step: One replaces 
time t everywhere by imaginary time r, t — > it, and restricts the integration range over r in the action 
to the interval [0, (3 = where T is the temperature of the system: 

Z = J VA^(x,T)V^(x,T)Vi>(x,T) e~ J'^ T J' d3xCE[A ^ ] . (2) 

The origin of this replacement is the realization that the partition function is defined as Z = tr p and that 
the density operator p = e~^ H for the grand canonical thermal equilibrium ensemble looks just like the 
time-evolution operator e tHt in the vacuum theory, with t replaced by i/3. The Euclidean Lagrangian 
density £e[A ip] arises from the normal QCD Lagrangian 

£qcd = ~\f; v F^ + (dp - i9Y A t) * ~ < 3 > 
with the non-Abelian gluon field strength tensor 

F%, = d„A a v - %A% + gfabcA^ , (4) 

by replacing dt = —id T as well as Aq = %A\ and jft = ij% (where = gipJu^fip is the color current 
vector of the quarks), and summing Lorentz indices over 1 through 4 (instead of through 1) with unit 

metric tensor. In addition, in order to preserve the invariance of Z[A, tp, ip] = tr p[A, ip, ip] under cyclic 
permutations of the field operators under the trace, the classical gluon fields in the path integral must 



be periodic in imaginary time, A^(x, r) = A^(x, t+(3), whereas the Grassmannian fermion fields obey 
antiperiodic boundary conditions, tp( x -> T ) = ~~ 4>{ x i T +0) et c Note that r is not really a time, it only 
plays a similar formal role in the path integral as real time does at zero temperature; a system in global 
thermal equilibrium is completely time independent. 

Similar to Eq. (2) one can write down a path integral for the thermal equilibrium ensemble average 
(O) = ti(pO) of an arbitrary observable 0[A, ip, ip] which depends on the quark and gluon fields: 

= J VA^jx, t) V$(x, t) IhHx, t) Q[A, j, e~ So drf^xC^M ^ 
J VA^(x, t) V$(x, t) V4>(x, t) e~ So dT I d 3 xC^[A,^] 

Here 0[A, ip, ip] is the classical observable, expressed through the classical fields A^(x, r), ip(x, r) and 
tp(x,r). 

So far this expression is exact. Since the QCD Lagrangian is bilinear in the quark fields, the path 
integral over the Grassmann field variables can be done analytically, resulting in an infinite-dimensional 
determinant over all space-time points: 



(0) = 



JVA^x,t) 0[A] det[i^ (d„ - ig%A*) - m] e~ So drfd 3 xC E [A] 



jVA^x.r) det[i 7 M - ig^A^ - m] e~ So drJd^c E [A] 



(6) 



Here £e[^4] = —\F^ V F£ V is the purely gluonic part of the euclidean QCD Lagrangian, and 0[A] arises 
from 0[A,ijj, ip] when doing the Gaussian integral over ip and -0. 

The words Lattice QCD stand for an algorithm to numerically evaluate this path integral, by 
discretizing space and time into N^N T space-time lattice points, evaluating the corresponding NgN T - 
dimensional fermion determinant, and integrating over the 4x8 gluon fields A a ^ from — oo to oo at each 
of the NgN T lattice points. The integrals are performed by Monte Carlo integration, using the Metropolis 
method of importance sampling. This method works well as long as the fermion determinant is positive. 
This is indeed the case for the expression given in Eq. (6) for the grand canonical ensemble at zero 
chemical potential. It describes a quark-gluon plasma with vanishing net baryon density, which is a good 
approximation for the Early Universe. Heavy-ion collisions, on the other hand, involve systems with 
non-zero net baryon number, brought into the collision by the colliding nuclei. This requires introduc- 
tion of a baryon chemical potential /j, b which enters into the fermion determinant in Eq. (6) with a factor 
i and leads to oscillations of the latter. The resulting "sign problem" has been a stumbling block for 
lattice QCD at finite net baryon density for almost 25 years, and only recently significant progress was 
made, resulting in first lattice QCD results for the hadronization phase transition in a quark-gluon plasma 
with nonzero baryon chemical potential [2] (although still limited to < 3T cr [3]). 

2.2 Color deconfinement and chiral symmetry restoration 

For our discussion two observables are of particular importance: the Polyakov loop operator 

(where A^ = A\^- is a 3x3 matrix and V stands for path ordering), and the scalar quark density 
-ip(x)ip(x). Due to translational invariance of the medium both have ^-independent thermal expecta- 
tion values which, however, show a strong temperature dependence. This is shown in Fig. 1. The 
argument of the exponential function in the Polyakov loop operator L is gauge-dependent although L 
itself is not (due to the trace and the periodicity condition A1(x,f3) = A%(x,0)); we can thus choose 
a gauge in which its r-dependence vanishes. It can then be interpreted as the interaction energy of 
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Fig. 1: Left: Polyakov loop expectation value (L) and its temperature derivative (Polyakov loop susceptibility \l) as a func- 
tion of the lattice coupling (3 = 6/g 2 which is monotonically related to the temperature T (larger (3 correspond to larger T). 
Right: The chiral condenstate an d the negative of its temperature derivative (chiral susceptibility Xm) as a function of 

temperature. (FromRef. [4].) 



an infinitely heavy quark at position x (whose euclidean color current density four-vector is given by 
j«( y ) = ig%5(y-x) (1, 0, 0, 0)) with the gluon field A£(y): 

where 

#int = ~L hlt = W d 3 yJ^(y)A^y) = J d 3 y JHy)A%(y) = ig^A%{x). (9) 

A vanishing thermal expectation value (L) of the Polyakov loop operator thus indicates infinite energy 
for a free quark, i.e. quark confinement. The left panel of Fig. 1 shows this to be the case at small 
temperatures. However, as the temperature increases, (L) increases rapidly to a nonzero value at high 
temperatures, with a relatively sharp peak of its derivative at a critical coupling f3 cr . This indicates that 
quark confinement is broken at the corresponding critical temperature T CT . 

In the right panel of Fig. 1 we see that at low temperatures the scalar quark density has a nonva- 
nishing expectation value ("chiral condensate") which evaporates above a critical coupling. Again the 
corresponding susceptibility shows a relatively sharp peak, at the same value (3 CV . In the absence of quark 
masses the QCD Lagrangian is chirally symmetric, i.e. invariant under separate flavor rotations of right- 
and left-handed quarks. Since the up and down quark masses in £qcd are very small, neglecting them 
is a good approximation. The nonvanishing chiral condensate at T = breaks this chiral symmetry and 
generates a dynamic mass of order 300 MeV for the quarks; the corresponding "constituent" masses in 
vacuum are thus about 300 MeV for the up and down quarks and about 450 MeV for the strange quark 
(whose bare mass in £qcd is already about 150 MeV). According to the right panel of Fig. 1 the dy- 
namically generated mass melts away at T cr , making the quarks light again above T cr : the approximate 
chiral symmetry of QCD is restored. 

Obviously deconfinement and chiral symmetry restoration happen at the same critical tempera- 
ture T cr . Figure 1 shows this for one specific, temperature-dependent value of the quark mass used 
in the lattice calculation (m q = 0.8 T). This value is unrealistically large, but calculations with realis- 
tic and temperature-independent masses are very costly and not yet available. Instead, one repeats the 
calculations for several unrealistically large masses and tries to extrapolate to zero mass. The perfect 
coincidence of the peaks in the chiral and Polyakov loop susceptibilities is seen for all quark masses [4] 
and thus expected to survive in the chiral limit. 



Both deconfinement and chiral symmetry restoration are phenomenologically important. Decon- 
finement leads to the liberation of a large number of gluons which can produce extra quark-antiquark 
pairs and drive the system towards chemical equilibrium among quarks, antiquarks and gluons. The melt- 
ing of the dynamical quark masses above T cr makes the quarks lighter and lowers the quark-antiquark 
pair production threshold. This is particularly important for strange quarks whose constituent quark 
mass is much higher than the critical temperature while its current mass is comparable to T cr . Above T cr 
thermal processes are therefore much more likely to equilibrate strange quark and antiquark abundances 
during the relatively short lifetime time of a heavy-ion collision. 

The dissolution of massive hadrons into almost massles quarks and gluons at T cr leads to a very 
rapid rise of the energy density near the deconfinement transition. This is shown in Fig. 2. For a massless 
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Fig. 2: Energy density in units of T 4 for QCD with two and three dynamical quark flavors [5, 6]. (The symbol e is here used 
for the energy density e.) The curves labelled "2 flavour" and "3 flavor" were calculated for two and three light quark flavors 
of mass = 0.4, respectively. "2+1 flavour" indicates a calculation for two light and one heavier strange quark flavor of 
5? = 1. In this case the ratio e/T 4 interpolates between two light flavors at T < T cr and three light flavors at T > 2T cr [6]. For 
2 and 2+1 flavors the critical temperature is T CI = 173 ± 15 MeV [6]. The right figure indicates, for the case of 3 light flavors, 
absolute values for the energy density e at several temperatures. 

gas of quarks and gluons the energy density is proportional to T 4 . The proportionality constant reflects 
the number of massless degrees of freedom, multiplied by ^ for bosons and by ||q for fermions. The 
arrow in the upper right of Fig. 2 indicates this constant evaluated for 2(helicity)x8(color)=16 gluon 
degrees of freedom plus 2(spin)x3(color)x3(flavor)x2(g+g)=36 massless quark and antiquark degrees 
of freedom. We see that for T < 4T cr the lattice data remain about 20% below this Stefan-Boltzmann 
limit. Near T cr the ratio e/T 4 drops rapidly by more than a factor 10. This is due to hadronization. 
The much heavier hadrons are exponentially suppressed below T cr , leading to a much smaller number of 
equivalent massless degrees of freedom. (Due to the somewhat too large quark masses in the simulations 
presented in Fig. 2 this effect is slightly exaggerated.) According to Fig. 2 the critical energy density for 
deconfinement is about 0.6-0.7 GeV/fm 3 . The ±15 MeV uncertainty in T cr induces a ±40% uncertainty 
in the critical energy density which could be as large as 1 GeV/fm or as small as 500MeV/fm 3 . I'll use 
e cr = 1 GeV/fm 3 as a ball-park number. It is easy to remember because it is the energy density in the 
center of a proton and corresponds to one proton mass in a cube of 1 fm along each side. 

The right panel of Fig. 2 also illustrates how painful the factor T 4 is: if we want to exceed the 
critical temperature by only 30% in order to reach the upper edge of the transition region, we already 
need an energy density e ~ 3.5 GeV/fm 3 , five times the critical value! And to reach 2T cr requires 
e ~ 23 GeV/fm 3 . The former number is approximately the value one was able to reach at the CERN 
SPS in Pb+Pb collisions at ^/snn = 17GeV. The latter value has been achieved at the Relativistic Heavy 
Ion Collider (RHIC) in Au+Au collisions at ^/snn = 130 GeV (see later). An estimate where the Large 
Hadron Collider (LHC) at CERN will take us in 2008 is indicated near the right edge of the Figure. The 
corresponding energy densities are of the order of 500 GeV/fm 3 ! 



2.3 The QCD phase diagram and how to probe it with heavy-ion collisions 

Figure 3 shows the phase diagram of strongly interacting matter in the temperature vs. baryon chemical 
potential plane (T, hb)- Cold nuclear matter, such as in the interior of, say, a Pb nucleus, sits a T = 
and fj,B ~ mN = 940 MeV. The short line emerging from this blob indicates the nuclear liquid-gas 
phase transition, with a critical endpoint at a temperature of about 7.5 MeV. At higher temperatures 
more and more hadron resonances are excited and we have a hadron resonance gas. From QCD lattice 
calculations we know the phase structure along the temperature axis, with the deconfinement transition 
from a hadron resonance gas to a quark-gluon plasma at T cr 170 MeV. Lattice QCD tells us that even 
for realistically small up and down quark masses the transition at \i% = is most likely not a sharp phase 
transition but a rapid crossover as shown in Fig. 2 [6]. Phenomenological models have long indicated [7] 
that at non-zero // b the QGP and hadron gas are separated by a critical line of roughly constant critical 
energy density e cr ~ 1 GeV/fm 3 . Improved QCD inspired models [8] and recent lattice calculations 
at moderate non-zero baryon chemical potential [2, 3] indicate that the transition becomes first order at 
non-zero [i b although the precise value where this happens still remains to be determined by performing 
calculations with more realistic smaller quark masses. At low temperatures and asymptotically large 
baryon densities quarks are also deconfmed, although not in a quark-gluon plasma state but rather in a 
color superconductor [8]. The superconducting state is separated from the QGP by a first order transition 
at a critical temperature estimated to be of order 30-50 MeV [8]. Whether it has a direct transition to 
normal nuclear matter as indicated in my sketch or if other phases (e.g. involving pion or kaon conden- 
sates) intervene is presently not known. There exists a rich spectrum of theoretical possibilities [8,9]. 




Fig. 3: The QCD phase diagram. 

Figure 3 also indicates a few typical trajectories through the phase diagram which we can follow 
in relativistic heavy-ion collisions. The collisions start at the cold nuclear matter point in the phase 
diagram, go through an early non-equilibrium stage which can not be mapped onto the phase diagram 
(indicated by dotted lines), and then reappear in the phase diagram after having thermalized at some high 
temperature. Unfortunately, one cannot use heavy-ion collisions to compress nuclear matter without 
producing a lot of entropy and therefore also heating it; hence it seems impossible to probe with them 
the color superconducting phase of strongly interacting matter. 

At the SIS (heavy-ion synchrotron) at the GSI facility in Darmstadt, with beam energies of order 
1-2 GeV/nucleon, the center of mass collision energy is not large enough to reach the deconfmed state 



of QCD. At the Brookhaven AGS (beam energies of about lOGeV/nucleon) we may have straddled the 
transition line. First clear indications that we had crossed into a new state of matter beyond the hadron 
resonance gas came from Pb+Pb collisions at 160GeV/nucleon beam energy [10]. Au+Au collisions at 
RHIC thermalize at an initial temperature To « 2 T cr (see later). By collective expansion the collision 
fireball cools down and passes through the hadronization phase transition from above. In this sense 
heavy-ion collisions do not probe the de-confinement, but rather the confinement phase transition, just as 
the early universe. As the center of mass collision energy increases, the colliding nuclei become more and 
more transparent, meaning that a decreasing fraction of the beam energy and of the incoming baryons 
get stopped in the center of mass system. The midrapidity collision fireball therefore contains fewer 
and fewer of the incoming net baryons, becoming more and more baryon-antibaryon symmetric (i.e. /j,b 
decreases). At RHIC the entropy per baryon ratio is between 200 and 300, depending on collision energy. 
At the LHC one expects S/A to go up to several thousand. While still far from the early universe ratio 
of 10 9 , this is, for all practical purposes, "baryon-free" (/xb = 0) QCD matter. 

3 THE DIFFERENT STAGES OF A HEAVY-ION COLLISION - THE LITTLE BANG 
3.1 The production of hard probes during the early collision stage 

Figure 4 summarizes the key stages of relativistic heavy-ion collisions: thermalization, expansion, and 
decoupling. In the very early collision stages, before the bulk of the quanta, which are created from the 
fraction of the beam energy lost in the collision, have time to rescatter, "hard" particles with either a large 
mass or large transverse momenta p± » 1 GeV/c are created. Their creation involves large momentum 
transfers Q 2 ~ p\ » 1 GeV 2 , therefore their production can be calculated in perturbative QCD, using 
factorization theorems, from the nuclear structure functions. According to the uncertainty relation hard 
particle production happens on a time scale Tf orm ~ 1/ ^/Q 2 ; for a 2 GeV particle this means Tf orm ~ 
O.lfm/c. 

At SPS energies and below, hard particles are essentially only produced in the primary collisions 
between the projectile and target nucleons; the bulk of the produced particles have transverse momenta 
below 2 GeV/c and are too soft to produce hard particles via secondary collisions. High-p^ j et s from the 
fragmentation of hard partons have practically relevant cross sections only at RHIC energies and above. 
Once produced, one can use them to probe the soft matter created by the bulk of soft particles [11]: 
In a central collision between two Pb or Au nuclei the nuclear reaction zone has a transverse diameter 
of about 12 fm, so a hard particle created near the edge and moving in the "wrong" direction (namely 
straight inward) needs 12fm/c before it emerges on the other side. During this time the soft matter 
thermalizes, expands, cools down and almost reaches decoupling, all of which is probed by the fast 
particle on its way to the other side. It does so by scattering off the evolving medium and losing energy 
which can be measured. The energy loss is proportional to the density of the medium times the scattering 
cross section between the probe and the medium constituents, integrated along the probe's trajectory. 

At SPS energies the "hardest" particles which can be produced and well identified by their dis- 
tinctive decay pattern are cc pairs. If produced in a proton-proton collision where either no hot soft 
matter is created or the charmed quarks leave the region in which soft particles are formed before that 
happens, these charmed quarks and antiquarks form either a bound charmonium state (J/tp, ip', or %, 
"hidden charm" production) or they find light quark partners to hadronize into "open charm" states (D 
and D mesons or charmed baryons). The corresponding branching ratios are well-known - the hidden 
charm states form a very small fraction (less than 1%) of all final states. If the cc pair is created in a 
heavy-ion collision, something similar as with the jets discussed above happens: the two heavy quarks 
have to travel through a dense medium of soft particles which interferes with their intention to hadronize 
and modifies their branching ratios into open and hidden charm states. In particular, if the soft medium 
thermalizes into a quark-gluon plasma, the color interaction between the c and c is Debye screened by 
the colored quarks and gluons in the plasma, thereby prohibiting their normal binding into one of the 
charmonium states. This should cause "J/Y> suppression" [12]. 
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Fig. 4: Stages of a relativistic heavy-ion collision and relevant theoretical concepts. 

Other probes of the early collision stage are direct photons [13], either real or virtual, in which 
case they materialize a lepton-antilepton (e + e~ or pairs, generally known as "dileptons". Such 

photons are produced by the electric charges in the medium (i.e. by quarks and antiquarks during the 
early collision stage). Their production cross section is proportional to the square of the Sommerfeld 
fine structure constant a = and thus small, but since they also reinteract only with this small elec- 
tromagnetic cross section their mean free path even in a very dense quark-gluon plasma is of the order 
of 50000fm and thus much larger than any conceivable heavy-ion fireball. In contrast to all hadronic 
probes, they thus escape from the collision zone without reinteraction and carry pristine information 
about the momentum distributions of their parent quarks and antiquarks into the detector. Real and vir- 
tual photons are emitted throughout the collision and expansion stage, and their measured spectrum thus 
integrates over the expansion history of the collision fireball. However, their production rate from a 
thermal system scales with a high power of the temperature (~ T 4 in a static heat bath but ~ T 6 once 
the effects of expansion are taken into account [13]) and is thus strongly biased towards the early col- 
lision stages. Unfortunately, the directly emitted photons and dileptons must be dug out from a huge 
background of indirect photons and and a large combinatorical background of uncorrected lepton pairs 



resulting from electromagnetic and weak decays of hadrons after hadronic freeze-out. This renders the 
measurement of these clean electromagnetic signals difficult. 

3.2 Thermalization and expansion 

The key difference between elementary particle and nucleus-nucleus collisions is that the quanta created 
in the primary collisions between the incoming nucleons can't right away escape into the surrounding 
vacuum, but rescatter of each other. In this way they create a form of dense, strongly interacting matter 
which, when it thermalizes quickly enough and at sufficiently large energy density, is a quark-gluon 
plasma. This is why with heavy-ion collisions we have a chance to recreate the matter in the very early 
universe whereas with high energy collisions between leptons or single hadrons we don't. 

The produced partons rescatter both elastically and inelastically. Both types of collisions lead 
to equipartitioning of the deposited energy, but only the inelastic collisions change the relative abun- 
dances of gluons, light and strange quarks. (To also change the abundance of the heavier charm quarks 
(m c ~ 1200 MeV) requires secondary collisions of sufficient energy which do not happen at the SPS but 
may play a role at RHIC and above.) From the phenomenology of pp collisions it is known [14, 15] that 
the produced hadron abundances are distributed statistically, but that strange hadrons are systematically 
suppressed (probably because strange quarks are not present in the initial state and their large constituent 
mass of ~ 450 MeV makes them hard to create from the vacuum as the produced partons hadronize). 
In a heavy-ion collision, if the reaction zone thermalizes at energy density > e cr such that gluons are 
deconfined and chiral symmetry is restored, strange quarks are much lighter (m s ~ 150 MeV) and can 
be relatively easily created by secondary collisions among the many gluons, leading to chemical equi- 
libration between light and strange quarks [16]. The observed strangeness suppression in pp collisions 
should thus be reduced or absent in relativistic heavy-ion collisions [17]. 

A thermalized system has thermal pressure which, when acting against the surrounding vacuum, 
leads to collective (hydrodynamic) expansion of the collision fireball. As a consequence, the fireball 
cools and its energy density decreases. When the latter reaches e cr ~ 1 GeV/fm 3 , the partons convert to 
hadrons. During this phase transition the entropy density drops steeply over a small temperature interval 
(similar to the energy density shown in Fig. 2). Since the total entropy cannot decrease this implies 
that the fireball volume must increase by a large factor while the temperature remains approximately 
constant. The growth of the fireball volume takes time, so the fireball ends up spending significant time 
near T cr . Furthermore, while the matter hadronizes its speed of sound c s = ^Jdpjde is small [6], causing 
inefficient acceleration so that the collective flow does not increase during this period. This may be 
visible in the direct photon spectrum: an inverse Laplace transform should show a particularly strong 
weight at the transition temperature, blueshifted by the prevailing radial flow as the system crosses T cr 
[18]. The softness of the equation of state near the phase transition should also manifest itself in the 
center-of-mass energy dependence of the collective flow, as discussed later. 

3.3 Hadronic freeze-out and post-freeze-out decays 

After hadronization of the fireball, the hadrons keep rescattering with each other for a while, contin- 
uing to build up expansion flow, until the matter becomes so dilute that the average distance between 
hadrons exceeds the range of the strong interactions. At this points all scattering stops and the hadrons 
decouple ("freeze out"). Actually, their abundances already freeze out earlier when the rates for inelastic 
processes, in which the hadrons change their identity, become too small to keep up with the expansion. 
Since the corresponding inelastic cross sections are only a small fraction of the total cross section, inelas- 
tic processes stop long before the elastic ones, leading to earlier freeze-out for the hadron abundances 
than for their momenta: chemical freeze-out precedes thermal or kinetic freeze-out. What I call "elastic" 
includes resonant processes such as tt+N — > A — > tt+N where two hadrons form a short-lived res- 
onance which subsequently decays back into the same particles (possible with different electric charge 



assignments). Such processes don't change the finally observed chemical composition, but they con- 
tribute to the thermalization of the momenta and have large "resonant" cross sections. Since most of the 
hadrons in a relativistic heavy-ion fireball are pions (since they are so light), resonances with pions are 
very efficient in keeping the system in thermal equilibrium (even after chemical equilibrium has been 
broken). For scattering among pions the p resonance plays a large role while kaons and (anti)baryons 
couple to the pion fluid via the K* , A and Y* resonances. Due to the particularly large A and Y* reso- 
nance cross sections, even at RHIC and LHC (where the net baryon density is small) baryon-antibaryon 
pairs play an important role as part of the "glue" that keeps the expanding pion fluid thermalized well 
below the chemical freeze-out point. 

At kinetic freeze-out all hadrons, including the then present unstable resonances, have an approxi- 
mately exponential transverse momentum spectrum reflecting the temperature of the fireball at that point, 
blueshifted by the average transverse collective flow. The unstable resonances decay, however, producing 
daughter particles with, on average, smaller transverse momenta. The experimentally measured spectra 
off stable hadrons cannot be understood without adding these decay products to the originally emitted 
spectra. Since most resonances decay by emitting a pion, this effect is particularly important for the pion 
spectrum which at low p± are completely dominated by decay products. This seriously affects the slope 
of their spectrum out to p± ~ 700 MeV [19, 20], making it steeper than the blueshifted thermal spectrum 
of the directly emitted hadrons. 

3.4 Theoretical tools 

The idea that the parton production process factorizes into a perturbatively calculable "hard" QCD cross 
section and a non-perturbative, experimentally determined nuclear parton structure function (or parton 
distribution function) can be used to describe the initial production of hard partons, with px > Po where 
Po > 1— 2GeV describes the lower applicability limit of this type of approach [21, 22]. The production 
of soft partons is usually non-perturbative and requires phenomenological models, such as string mod- 
els [23], for their description. At very high collision energies particle production at midrapidity (i.e. 
particles with small longitudinal momenta in the center-of-momentum frame) probes the nuclear struc- 
ture functions at small x where x is the fraction of the beam momentum carried by the partons whose 
collision produces the secondary partons. At small x the gluon distribution function becomes very big 
and gluons begin to fill the transverse area of the colliding nuclei densely [24, 25]. The gluons start to 
recombine, leading to gluon saturation, and low p^ gluon modes are occupied by a macroscopic number 
of gluons ~ l/a s where the effective strong coupling a s is small due to the high density of gluons [26]. 
As a result, the initial gluons can be effectively described by a classical gluon field [26, 27] in which the 
coupling is weak but nonlinear density effects are important (this state has become known as the "color 
glass condensate" [28]). The production of soft secondary gluons can be understood as the liberation 
of these gluons by breaking the coherence of their multiparticle wavef unction [27, 29]. The collision 
energies where these modern ideas become applicable lie probably beyond the RHIC range, but this is 
an exciting and active field of ongoing research which may come to fruition at the LHC. 

After the initial parton production or "liberation" process one must describe the rescattering and 
thermalization of the produced quanta. For not too dense systems this can be done with classical kinetic 
transport theory (relativistic Boltzmann equation), also known as parton cascade [30, 31]. In the early 
collision stages at RHIC and LHC energies the densities are probably too high for this to remain a 
reliable approach, and one must switch to quantum transport theory [32]. This formalism has not yet been 
developed very well for practical applications, and much interesting work is presently going on in this 
direction (e.g. [33, 34]). None of these approaches, which usually invoke perturbative QCD arguments 
to describe the microscopic scattering processes, has so far been able to produce parton thermalization 
time scales at RHIC which are shorter than about 5fm/c [33]. This is too long for heavy-ion collisions 
whose typical expansion time scale (Hubble time) is only a couple of fm/c. As I will show later there is 
strong phenomenological evidence that thermalization must happen much more quickly. This presents 



an interesting challenge for theory and indicates the importance of strong non-perturbative effects in the 
early collision stages and the QGP. 

Once local thermal equilibrium has been reached, the further evolution can be described hydrody- 
namically. The simplest version of such an approach is ideal fluid dynamics which will be discussed in 
more detail later. The fireball can be described as an ideal fluid if the microscopic scattering time scale 
is much shorter than any macroscopic time scale associated with the fireball evolution. If this is not the 
case, one should include non-ideal effects such as shear, diffusion and heat conduction. This requires the 
calculation of the corresponding transport coefficients in a partonic system close to thermal equilibrium. 
Much work is going on in this direction (see e.g. [35, 36, 37]), but systematic approaches based on 
perturbatively resummed thermal QCD [35, 36, 37] give phenomenologically unacceptable large values 
for these transport coefficients, again indicating the importance of non-perturbative effects. 

The hydrodynamic equations require knowledge of the equation of state, i.e. the relationship 
p(e,riB) between the pressure, energy density and baryon density. For small ub this is known from 
lattice QCD; for large rig one extrapolates the lattice results with phenomenological models (see [9] for 
a recent review). Hydrodynamics is the ideal language for relating observed collective flow phenomena 
to the equation of state. Since it only requires the equation of state but no detailed knowledge of the 
microscopic collision dynamics, it allows for an easy description of the hadronization phase transition 
without any need for a microscopic understanding of how hadrons form from quarks and gluons. Of 
course, the underlying assumption is that all these microscopic processes happen so fast that the system 
never strays appreciably from a local thermal equilibrium. 

After hadronization the system continues to expand and dilute until the average distance between 
hadrons becomes larger than the range of the strong interaction. At this point the hadrons decouple, i.e. 
their momenta stop changing until they are recorded by the detector. One can implement this decoupling 
or "kinetic freeze-out" in different ways, either by truncating the hydrodynamic phase abruptly with the 
Cooper-Frye algorithm (see below) or by switching back from hydrodynamics to a (this time hadronic) 
cascade [38, 39, 40] in which decoupling happens automatically and selfconsistently. In the Cooper-Frye 
algorithm one also must take into account that in this approach all kinds of hadrons decouple simulta- 
neously, and that unstable hadron resonances decay subsequently by strong (and in some situations also 
by weak) interactions before their stable daughter products reach the detector. To compute the measured 
spectra one must therefore fold the initially emitted Cooper-Frye hadron spectra with these decays [19]. 
Some studies, such as two-particle momentum correlations, also require the computation of long-range 
final state interaction effects, such as the Coulomb repulsion/attraction between charged hadrons which 
continues long after their strong interactions with each other have ceased. 

3.5 Strategies for reconstructing the Little Bang 

The bulk (over 99%) of the particles produced in heavy-ion collisions are hadrons. These are strongly 
interacting particles which cannot decouple from the fireball before the system is so dilute that strong 
interactions cease. Their observed momenta thus provide a snapshot of the kinetic decoupling stage 
("thermal freeze-out"). In this sense hadrons are the Little Bang analogue of the cosmic microwave 
background in the Big Bang. As discussed in Sec. 3.3, hadron abundances freeze-out earlier. As we will 
see later, this "chemical freeze-out" happens right after hadronization, i.e. the finally observed hadron 
abundances are generated by the hadronization process itself. This primordial hadrosynthesis is the Little 
Bang analogue of primordial nucleosynthesis in the Big Bang. As the hadrons decouple, they carry not 
only thermal information about the prevalent temperature at chemical resp. thermal freeze-out, but in 
the momentum spectra this information is folded with (i.e. blueshifted by) the collective expansion flow, 
just as the temperature of the cosmic microwave radiation is redshifted by the cosmological expansion. 
The hydrodynamic expansion flow is the Little Bang analogue of the cosmic Hubble expansion in the 
Big Bang. (Of course, the origin of the expansion is entirely different in the two cases: in heavy-ion 
collisions it is generated hydrodynamically by pressure gradients whereas in the Big Bang it reflects an 



initial condition, modified over billions of years by the effects of the gravitational interaction. But the 
velocity profiles turn out to be surprisingly similar, as we will see!) 

We reconstruct the Little Bang from these "late" hadronic observables very much like we recon- 
structed the Big Bang from the 3 pillars of cosmology, Hubble expansion, CMB, and primordial nuclear 
abundances. The hadronic observables are abundant and can be measured with high statistical accuracy. 
As I will show, their theoretical analysis allows to separate the thermal from the collective motion. The 
collective flow provides a memory of the pressure and other thermodynamic conditions during the ear- 
lier collision stages. In fact, we will see that certain anisotropics in the flow patterns seen in non-central 
heavy-ion collisions ("elliptic flow") provide a unique window into the very early collision stages and 
are no longer changed after about 5 fm/c after initial impact. As you will see, a (in my opinion) water- 
tight proof for thermalization during this early stage, at a time T t h C rm < 1 fm/c and at prevalent energy 
densities which exceed the critical value e cr for deconfmement by at least an order of magnitude, can 
be based on the accurately measured elliptic flow of the final state hadrons. (This is a bit analogous to 
the indelible imprint that cosmic inflation has left on the density fluctuations in the universe, which can 
be accurately measured through the anisotropy of the cosmic microwave radiation which only decoupled 
300 000 years later.) 

The reconstruction of the global space-time evolution of the reaction zone from the finally ob- 
served soft hadrons is the cornerstone of the program. The reconstructed dynamical picture will be the 
basis on which other rarer observables, in particular the "deep" or "hard" probes, will be interpreted. For 
example, jet quenching and heavy quarkonium suppression cannot be quantitatively interpreted without 
knowledge of the fireball density and its space-time evolution, and direct photon and dilepton spectra 
cannot be properly understood without a relatively accurate idea about the transverse flow patterns at the 
time of photon emission. Still, these probes will in the end be the only direct access we will ever have 
to the temperature and energy density at the beginning of the expansion when the system was a quark- 
gluon plasma, and although they cannot be fully interpreted and exploited without the later emitted soft 
particles they are still an indispensable part of the picture. This illustrates the network-like interdepen- 
dence between soft and hard observables in their role for elucidating heavy-ion collision and quark-gluon 
plasma dynamics. 

4 THE LITTLE BANG - COLLECTIVE EXPLOSION OF A THERMALIZED SYSTEM 

An unavoidable consequence of quark-gluon plasma formation in heavy-ion collisions is collective flow. 
Since a quark-gluon plasma is by definition an (approximately) thermalized system of quarks and gluons, 
it has thermal pressure, and the pressure gradients with respect to the surrounding vacuum cause the 
quark-gluon plasma to explode. Absence of collective flow would indicate absence of pressure and 
imply absence of a hot thermalized system and, a fortiori, of a quark-gluon plasma. 

In this chapter we'll therefore do two things: (i) learn how to analyze the measured particle spectra 
for the presence of collective flow and how to separate it from random thermal motion, and (ii) compute 
the expected collective flow patterns from reasonable initial conditions using a hydrodynamic model. By 
comparing with experiments we can fine-tune the initial conditions, learn about the equation of state of 
the hot expanding matter and, most importantly, find out about when and at which energy densities the 
thermal pressure builds up and begins to drive the collective expansion. 

4.1 Radial flow 

4.1.1 Flow defined 

Consider a nuclear fireball undergoing collective expansion. Collective flow is defined by the following 
operational procedure: at any space-time point x in the fireball, we consider an infinitesimal volume el- 
ement centered at that point and add up all the 4-momenta of the quanta in it. The total 3-momentum P 
obtained in this way, divided by the associated total energy P°, defines the average "flow" velocity v(x) 



of the matter at point x through the relation P/P° = v. Collective flow thus describes a correlation be- 
tween the average momentum of the particles with their space-time position, a so-called x-p-correlation. 
With v{x) we can associate a normalized 4-velocity = 7(1, v) where j(x) = 1/^/1— v 2 (x) is the 
corresponding Lorentz dilation factor and u-u = u^u^ = 1 (in units where the speed of light c = 1). 

I separate the flow velocity v{x) into its components along the beam direction ("longitudinal 
flow" v l) and in the plane perpendicular to the beam ("transverse plane") which I call "transverse flow" 
v±. The magnitude v± may depend on the azimuthal angle around the beam direction, i.e. on the 
angle between vj_ and the impact parameter b of the collision. In this case we call the transverse flow 
"anisotropic". Its azimuthal average we call radial flow. 



4.1.2 Local thermodynamic equilibrium 

If the fireball is in local thermodynamic equilibrium, we can not only define a local flow 4-velocity 
u^{x), but also a local temperature T(x) and, for each particle species i, a chemical potential 
which controls its particle density at point x. In this case the phase-space distribution of particles of type 
i is given by the Lorentz covariant local equilibrium distribution 
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Here gi is a spin-isospin-color-flavor-etc. degeneracy factor which counts all particles with the same 
mass mi and chemical potential fii. The factor p ■ u(x) in the exponent is the energy of the particle in the 
local rest frame (local heat bath frame), boosted to the observer frame by the flow 4-velocity of 
the fluid cell at point x (p ■ u^p° = E when — > (1, 0)). The ±1 in the denominator accounts for the 
proper quantum statistics of particle species i [upper (lower) sign for fermions (bosons)]. The Boltzmann 
approximation corresponds to keeping only the first term in the sum over n in the last expression. In our 
applications this is an excellent analytical approximation for all hadrons except for the pion. Pions and 
quarks and gluons are too light, mi < T, and one must use the proper quantum statistical distributions. 



4.1.3 Rapidity coordinates 

At relativistic energies it is convenient to parametrize the longitudinal flow velocities and momenta in 
terms of rapidities (for any velocity v the associated rapidity is 77 = \ In j±| or v = tanh rj): 



7?t = — In , y = — In = — in . (11) 
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As «l — > 1, rji,—> 00. Rapidities have the advantage over longitudinal velocities that they are addi- 
tive under longitudinal boosts: a fluid cell with flow rapidity 7?l in a given inertial frame has rapidity 
7]' L = rjL + Ar] in another inertial frame which moves relative to the first frame with rapidity A77 in the 
—z direction. The flow 4- vector (u°, u) = (u°, u±,ul) is then parametrized as 

= 7-L( cosnr /L, v x , v y , smhiji) with 7_l = = , (12) 

where = (v x , v y ) is the transverse flow velocity. Similarly the 4-momentum p^ = (E,p±,pi,) is 
written as 

p^ = (m±coshy, p x , p y , m^sinhy) with m± = \J m 2 +p 2 L = ^Jm 2 +p 2 +p 2 (13) 



which obviously satisfies the mass-shell constraint p 2 =p fl p fl = m 2 . [We use the notations p±=pt 
interchangeably for the transverse momentum, and similarly tbt = m± for the "transverse mass" - in 



the literature one often also finds the notations p t and m t .] The scalar product p ■ u(x) in the exponent 
of the Boltzmann factor then becomes 

p-u = j±(m ± cosh(y-7] L )-v ± -p ± y (14) 

4.1.4 Longitudinal boost-invariance and Bjorken scaling 

So far no approximations have been made as long as we allow v± and 7?l to be arbitrary functions of 
= (t, r±, z) where r± = (x, y) denotes the transverse coordinates. Things simplify, however, enor- 
mously if one assumes longitudinal boost-invariance. Bjorken [41] argued that at asymptotically high 
energies the physics of secondary particle production should be independent of the longitudinal reference 
frame. This condition can be easily expressed if one puts the nuclear collision at longitudinal position 
z = and introduces, instead of z and t, the following "space-time rapidity" and "longitudinal proper 
time" coordinates to describe the forward light cone emanating from the collision point: 

1 t+Z 1 1 + 1 rr; k 

V=2 ln — z = 2 ln TI> ^ = 7^2. (15) 

In these coordinates x M = (t, r±,z) reads x M = (rcosh?7,r^,rsmh77), and the space-time integration 
measure takes the form d A x = r dr di] d 2 r±. In a i-z-diagram, lines of constant space-time rapidity 
are rays through the origin with slope 1 / tanh rj. Longitudinal boost-invariance of particle production 
then implies that the initial conditions for local observables (such as particle and energy densities) are 
only functions of r and r±, but independent of r\ [41]. Furthermore, the boost-invariance of these 
initial conditions is preserved in longitudinal proper time if the system expands collectively along the 
longitudinal direction with a very specific "scaling" velocity profile = f [41]- Inserting this into the 
definition (11) for the longitudinal fluid rapidity t/l and comparing with (15) we see that longitudinal 
boost-invariance implies the identity 

7?l = rj (Bjorken scaling) (16) 

of the longitudinal fluid rapidity with the space-time rapidity. Since ^(x) characterizes the average 
longitudinal momentum of the produced particles at point x whereas rj characterizes the coordinate x 
itself, this identity implies a very strong correlation between the average longitudinal momentum and the 
longitudinal position. 

The Bjorken scaling approximation is expected to be good at high energies and not too close to the 
beam and target rapidities, i.e. in safe distance from the longitudinal kinematic limits. We will use it to 
describe particle production near midrapidity (y m in the cm. frame). As a result of Bjorken scaling, 
the Boltzmann exponent reduces to 

p-u(x) = 7±{r±,T)(m ± cosh(y-7]) - p ± • w±(rj_,r)). (17) 

Due to the additivity of rapidities, this is manifestly invariant under longitudinal boosts since it only 
involves rapidity differences. The transverse flow velocity, as well as the temperature T and chemical 
potential p,i, are independent of r\ and depend only on the transverse position r_|_ and longitudinal proper 
time r. 

4.1.5 The Cooper-Frye formula 

Suppose we want to count the total number of particles of species i after produced in the collision. Since 
this number does not depend on the reference frame of the observer, we must be able to express it in a 
Lorentz-invariant way. We define a three-dimensional hypersurface S(x) in 4-dimensional space-time 
along which we perform the counting. The simplest case would be a measurement in all space at a fixed 



global time t. In a space-time diagram this would correspond to a horizontal line at fixed t. But this is not 
how a real detector works. In an ideal detector, we surround the collision region hermetically by detector 
elements which sit stationary at a fixed distance from the collision point, and we wait until the particles 
pass through these detector elements (which particles with different velocities will do at different times). 
Assuming, for example, a spherical detector of radius R, the corresponding detection hypersurface would 
in a (t, r) space-time diagram be represented by a vertical line at fixed r = R, extending from t = — oo 
to t = oo. In this case the hypersurface is a 2-dimensional sphere extending over all time, i.e. it is again 
3-dimensional. 

You see that different choices for the 3-dimensional hypersurface S are possible as long as it 
completely closes off the future light cone emerging from the collision point. We count particles crossing 
the surface by subdividing it into infinitesimal elements d 3 a, defining an outward-pointing 4-vector 
d 3 (Tfj,(x) perpendicular to S(x) at point x with the magnitude d 3 a, computing the scalar product of 
the 4-vector j^(x) describing the current of particles % through point x, and summing over all such 
infinitesimal hypersurface elements: 

N t = J^d 3 o„{x)tf{x) = J^a,{x) J ^p"/i(*,p)) . (18) 

The particle number current density (x) is given in terms of the Lorentz-invariant phase-space distri- 
bution (giving the probability of finding a particle with momentum p at point x) by multiplying it with the 
velocity ^ and integrating over all momenta with measure ^ = ^h)* wriere h = 2nh = 2tt is Planck's 
quantum of action and we are using units where Ti = 1. 

Dividing by the Lorentz-invariant momentum-space measure ^ (which in rapidity coordinates 
reads dyp±dp± dip p = dy m±dm± dip p where ip p is the azimuthal angle of p±) we obtain the invariant 
momentum distribution for particle species i: 

E ^i = ^ = d Jfl = J- [ p .d*o (x)fAx,p). (19) 
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This is the Cooper-Frye formula [42]. One can show that two different surfaces Si and £2 give the 
same particle number Ni if between Si and £2 the distribution function fi(x,p) evolves via a Boltzmann 
equation with a collision kernel which preserves the number of particles i, and that we obtain the same 
momentum spectrum if and only if fi(x,p) evolves from Si to S2 by free-streaming, i.e. if fa is a 
solution of the collisionless Boltzmann equation. To compute the measured momentum spectrum we 
can therefore replace the surface S corresponding to the detector by shrinking it to the smallest and 
earliest surface that still encloses all scattering processes. We call this the "surface of last scattering" or 
"freeze-out surface" Sf. 

Any 3-dimensional space can be parametrized by a set of three locally orthogonal coordinates 
u, v, w. The points on the surface then have coordinates S M (u, v, w), p = 0,l, 2, 3. The normal 
vector d 3 a^(x) is computed through the formula 

, do v da x da? 

d a^ = -e^xp ^— ^— du dv dw (20) 
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where e^xp is the completely antisymmetric Levi-Civita symbol in four dimensions, with e 0123 = 
— eoi23 = 1- If we assume longitudinal boost-invariance, the freeze-out surface Sf can be character- 
ized by a longitudinal proper time Tf(r±) such that 

^i(r±,v) = {tf,Xf,Vf,Z{) = (r f (r ± ) cosh 77, r ± ,r f (r ± )sinh 77). (21) 

In this case 

d 3 ap= (coshrj,- — ,- — ,-smhr])T i (r ± )d 2 r ± dri, (22) 



and the integration measure in the Cooper-Frye formula (19) becomes (with (d x , d y ) = VjJ 



p ■ d 3 a(x) = (m± cosh(y— 77) — p± • Vj_Tf (r±)j Tf(r±) d 2 r± drj. 



(23) 



4.1.6 Thermal spectra from an exploding source 

The local thermal equilibrium form (10) for the phase-space distribution is a good approximation as long 
as the particles rescatter intensely, i.e. as long as the relaxation time for returning to equilibrium is shorter 
than any macroscopic time scale which changes the parameters T(x), fa(x) and u^(x). Use of Eq. (19) 
to compute the measured momentum spectrum requires knowledge of the phase-space distribution on 
the surface of last scattering, i.e. when all scattering has ceased. How can we relate the two? If the 
system expands very fast, its density decreases rapidly and the mean free path of the particles growth 
quickly. The transition from strong coupling to free-steaming thus happens in a short time interval. 
During this short time it is unlikely that the phase-space distribution undergoes qualitative changes, and 
we may approximate fi(x,p) on the last scattering surface by its thermal equilibrium form that it still 
had just a little earlier. This results in the famous Cooper-Frye freeze-out algorithm: determine the 
surface where local equilibrium begins to break down and use the local equilibrium distribution (10), 
with parameters T(x), fa(x) and u^(x) determined on this surface, in Eq. (19) to compute the final 
momentum distribution. The result is [43, 44] (note that T{, T, fa, v± and j± are all functions of r±) 
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where 0±(r±) = m± ^(rt) ' ^ ue to ^ e assumed boost-invariance, which allowed for an easy integra- 
tion over 77, this result is independent of rapidity y. 

For central collisions (6 = 0) the fireball is azimuthally symmetric, so v± = v±e r and (in polar 
coordinates r± = (r±, ip s )) the functions Tf, T, pi, v±, 7j_ are all independent of ip s . This also implies 
Vj_Tf = (dTi/dr±)e r± . The only angular dependences then arise from the scalar products v±-p± and 
Px-ViTf which both involve the relative angle between p± and r±, i.e. cos(v? s — ip p ). The azimuthal 
integral can thus be done analytically, producing another set of modified Bessel functions: 
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x(m ± K 1 (n(3±)I (na±) - p±-7^K (n(3 ± ) h(na±fj , (25) 

where a± = 7x ^ p± = a±(r±). For all hadrons except pions this can be used in the Boltzmann approx- 
imation, by keeping only the term n = 1. The factor Tf(r±) e^ i ^ r± ^ T( - r±S) =rii(r±) can be interpreted 
as the (unnormalized) radial density profile of the particles i. Introducing the radial flow rapdity p via 
v± = tanh p, which allows to write f3± = m ± c ° sh P and a± = p±s j, nhp , we then obtain the following 
"flow spectrum" [44,45]: 
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dy m±dm± 



^J r±drj.n i {r ± )^n A K 1 ( ^ )l ( ^ j (26) 

9t{ /m± cosh p{r±)\ fp±smh.p(r±) 
~ P± d^ K °\ T(r ± ) M T (r ± ) 



This formula is useful because it allows to easily perform systematic studies of the influence of the radial 
profiles of temperature, density and transverse flow on the transverse momentum spectrum, in order to 
better understand which features of a real dynamical calculation of these profiles control the shape of the 
observed spectra. There are many such case studies documented in the literature (see, e.g., [45, 46]), and 
I will here discuss only the most important and generic characteristics. 

4.1. 7 How radial flow affects single-particle transverse momentum spectra 

Since for all hadrons m±/T > 1, the modified Bessel functions K u can be approximated by exponentials 
~ e~ m± cosh P/ T _ The temperature on the freeze-out hypersurface is approximately constant [47, 48] 
since freeze-out is controlled by the mean free path which is inversely proportional to the density, which 
itself is a steep function of temperature [49]. Nevertheless, the flow spectra are characteristically curved, 
due to two effects: the influence of the l u Bessel functions at low p± and the integration over the radial 
flow profile p(r±). 

Let us first see what kind of flow profiles we should consider. At r± = the radial flow velocity 
must vanish by symmetry; as you follow the freeze-out surface out to larger r±, v± typically rises 
linearly with r± [47, 48, 39]. As shown in the right panel of Fig. 5, it eventually reaches a maximum 
value and drops again to zero since the dilute tail of the initial density distribution freezes out early 
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Fig. 5: Left panels: Hydrodynamical freeze-out surface Tf(r±) for non-central Pb+Pb collisions at & = 8fm (b points in x 
direction) with two different initial energy densities, one corresponding to SPS energies (left), the other corresponding to LHC 
energies (middle) [50]. Right panel: Radial flow rapidity profile p(r±) = yT(r±) for central Au+Au collisions at RHIC, from 
hydrodynamic calculations employing three different equations of state (for details see Ref. [39] from where this figure was 
taken). 

before radial flow develops. The left two panels of Fig. 5 show typical freeze-out surfaces Tf (rj_) from 
hydrodynamic calculations (see later) for SPS (left) and LHC energies (middle) [50]. One sees that 
at SPS energies the freeze-out surface moves from the edge inward since the fireball matter cools and 
freezes out faster than the developing radial flow can push it out. At LHC energies things begin similarly, 
but then the much stronger radial flow generated by the much higher internal pressure makes the fireball 
grow considerably before suddenly freezing out after about 13 fm/c. The freeze-out surface for RHIC 
energies lies in between [51], with transverse flow almost exactly balancing freeze-out and leading to an 
almost constant freeze-out radius (also seen in the right panel of Fig. 5) until after about 12 fm/c it rather 
suddenly shrinks to zero, again indicating sudden bulk freeze-out. 

So how does the radial flow affect the spectra? Looking at Eq. (refflowspec2) we see that in the 
absence of flow (p = 0) the second term vanishes (since li(0) = 0) and the first term simply becomes (for 
constant decoupling temperature) 



dy mj_dmj 



(27) 



This depends only m± (rather than on both m± and p±_ = y mj_— m?) and, except for the fact that m± 
cannot be less than the rest mass, the spectrum is therefore identical for all hadrons! This is known as 
"m±-scaling": In a static fireball all hadron spectra follow the same exponential dNi/(dym±dm±) ~ 
m 1 / 2 e~ m± / T , and the fireball temperature can be immediately extracted from their slope. 

Collective flow breaks mj_-scaling. This is shown in Fig. 6. The breaking of m^-scaling does 




Fig. 6: Flow spectra for various hadrons as a function of m±—mo where mo is their rest mass. The calculation assumes an 
infinitesimally thin shell of temperature T — 150 MeV expanding with v± — 0.4 (left) and tii = 0.9 (right). The curve labelled 
'V + (all)" includes pions from resonance decays in addition to the thermally emitted pions. 

not occur at large m± where again all spectra approach exponentials with identical slope. However, 
the inverse slope of these spectra now reflects a blueshifted freezeout temperature: Approximating 
p±~ m±^> mi and combining the Ko, Ki and Ii, Iq Bessel functions gives dNi/(dym±dm±) ~ 

exp[— ^(coshp— sinhp)] = exp(— mj_/T s i opc ) with the inverse slope parameter T s i ope = Ty / :|^-j- . 
This is a standard relativistic blueshift factor reflecting the boost of the thermal radiation towards the 
detector with radial flow velocity v±. 

Transverse flow breaks m± -scaling at small transverse momenta p± < mi where momenta and 
velocities can be added non-relativistically. The exact form of the of this breaking of m^-scaling depends 
on the density and velocity profiles. It is most extreme for a thin shell expanding with fixed velocity 
("blast wave"), shown in Fig. 6, in which case for sufficiently large hadron mass and flow velocity the 
spectrum develops a "blast wave peak" at nonzero transverse momentum [52]. More realistic calculations 
take into account the integration over a velocity profile in which the hole at low p± in the blast wave 
spectrum is filled in by contributions with smaller radial boost velocities from the fireball interior [45]. 
For a Gaussian density profile combined with a non-relativistic linear transverse velocity profile the 
integrated spectrum can be calculated analytically [53], and one finds again an exponential m± -spectrum, 
but this time with inverse slope Tj s i ope =Tf + ^mi(v±) 2 . We summarize these two important limits: 

• Non — relativistic, p± <C mi : T jS i ope Tf + -mi(v±) 2 (28) 



Relativistic, p± 3> m, : T s i ox > c Tf J — — — (for all mi) (29) 

V 1_V -L 



The slope systematics for hadrons with different masses in the low-p^ region thus allow to separate 
thermal from collective flow motion. Note that Eq. (28) can never be applied to pions since in the region 



p± -C m n their slope is affected by Bose statistics and by the contamination from resonance decays (see 
Fig. 6) neither of which is accounted for by Eq. (28). 



4.1.8 Extracting the freeze-out temperature and flow from measured transverse momentum spectra 

Equations (28) and (29) show that, as long as mi(v±) > 2Tf, the spectra are steeper at high p± and bend 
over becoming flatter at low p±. It is therefore difficult to characterize them by a single slope, especially 
when the detector measures different hadrons in different pj_ -windows, or when two different experi- 
ments measure the same hadron in different p^-windows. To extract the flow velocity using Eq. (28) 
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Fig. 7: Positively and negatively charged hadron spectra from central Pb+Pb collisions at 40, 80 and 160 A GeV beam energy 
at the SPS (left to right), measured by the NA49 Collaboration [54]. Also shown are 2-parameter fits with Eq. (30), assuming 
a sharp transverse flow velocity /3t (i.e. w(p) = S(p— tanh -1 /3t)) and constant rt . The resulting fit values for T and (3t 
are given in the figures. Positive and negativ hadrons were fitted independently; dashed lines indicate hadrons which were not 
included in the fit. The dashed pion curve does not include resonance decay contributions and Bose statistics. 



requires measuring all hadrons in a common interval of nonrelativistic transverse kinetic energy satis- 
fying m±—mi < rrii. Since such a procedure throws away information outside the common window, it 
is not very efficient. A much preferred method is to use the entire experimentally available information 



on the spectra by performing a simultaneous fit to all hadrons over all m± using Eq. (26). Assuming a 
constant freeze-out temperature, a common shape for the density profiles rii(r±), and a specific shape 
for T{(r±), we can rewrite Eq. (26) as an integral over flow rapidities [46], 



dy m±dmj 



Ni / dpw(p) 



n 

m±Kx\- 



mj_coshp\ /p±s'mhp 



T 



dr t /m_LCOshp 

- p± d^H—T— ]l1 



T 

p± sinh p 



(30) 



with the flow rapidity distribution 



Niw(p) = - J or ± (p)n i (r ± (p)) 



7T- 



dr± 
dp 



(3D 



where r±(p) is the inverse of the velocity profile p(r±). We see that for the shape of the spectrum the 
density and velocity profiles rii(r±) and vj_(r±) have no independent relevance - only the flow rapidity 
distribution w(p) matters. If we fix a reasonable shape for w(p) (see [46]), leaving its mean (p) free, we 
end up with a 2-parameter fit in terms of the freeze-out temperature T = Tf and average transverse flow 
velocity (v±) = tanh -1 (p). 

Figure 7 shows such a fit to hadronic m± -spectra measured by the NA49 Collaboration [54] in 
Pb+Pb collisions at the SPS at three different beam energies. Even disregarding the fitted lines the flow- 
typical hierarchy of slopes in the low-m^ region is obvious. The fit was done with a thin-shell model 
with fixed radial flow velocity. As discussed above this tends to exaggerate the bending of the spectra at 
low especially for the heavier particles 3, Q and the deuteron. This tendency can be seen in Fig. 7. 
Otherwise the quality of the fit, its ability to describe all spectra simultaneously, and the consistency of 
the fit parameters extracted from positively and negatively charged hadrons is impressive. 
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Fig. 8: Pion, kaon and antipro- 
ton spectra from 200 A GeV cen- 
tral Au+Au (left) and minimum 
bias p+p collisions (right), mea- 
sured by the STAR experiment 
[55]. Note the similar slopes for 
kaons and antiprotons in p+p col- 
lisions and their dramatically dif- 
ferent slopes at low transverse ki- 
netic energy in central Au+Au 
collisions. 
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The flattening of the spectra at low transverse kinetic energy m±—mo by transverse collective 
flow is even more dramatic at RHIC. Figure 8 shows a direct comparison of the negative pion, kaon and 
antiproton spectra in central Au+Au and minimum bias proton-proton collisions, at the same center of 
mass energy. Clearly, in p+p collisions the kaon and antiproton spectra have the same slope, indicating 
the absence of transverse collective flow. That the pion spectra are steeper than both kaons and antipro- 
tons can be attributed to the contribution of resonance decay pions which accumulate at low transverse 
momenta. However, the pion spectra in Au+Au are obviously flatter in Au+Au than in p+p, and this is 
even more true for kaons and antiprotons, with a large difference in slope between those last two. Even 
without a quantitative fit this is a clear manifestation of strong radial flow. 

A two-parameter flow fit on RHIC data from 200 A,GeV Au+Au collisions was performed by J. 
Burward-Hoy [56], see Fig. 9. Here a box-profile for the transverse density rii(r±) and a linear transverse 
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Fig. 9: Kinetic freeze-out temperature Tf and 
transverse flow velocity /3r at the fireball edge, 
extracted from a simultaneous fit of a flow spec- 
trum parametrization to n , K ± ,p and p spec- 
tra from 200 A GeV Au+Au collisions over the 
entire range of collision centralities.[56] More 
peripheral collisions (small numbers Af part of 
participating nucleons) are seen to decouple 
earlier, at higher freeze-out temperature and 
with less transverse flow. 
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velocity profile were used where the number /Jp given in Fig. 9 is the surface velocity (the average 
transverse velocity is (v±) = §/?t)- Again in the most central collisions freeze-out temperatures of about 
120 MeV and average radial flow velocities (v±) , ~ 0.45 are found. One should note, however, that in 
this case the pions were included in the fit, without correction for resonance decay feeddown. Resonance 
decay contributions steepen the pion spectrum, biasing the fit towards lower temperatures and/or flow 
velocities. The bands in the Figure attempt to provide an estimate of the systematic uncertainty associated 
with details of the model assumptions on the density and velocity profile and of the p± range of the pion 
spectrum used in the fit [56]. 

An interesting aspect of Fig. 9 is the centrality dependence of the fit parameters: The fits were per- 
formed over a very wide range of collision centralities, and one observes that more peripheral collisions 
tend to develop less radial flow and freeze out at a higher temperature. This is consistent with the ex- 
pectation that the collision fireballs in peripheral collisions don't live as long and have less time to build 
up radial flow. However, even though their average flow is smaller, the flow gradients are bigger (due to 
the smaller size), leading to decoupling at higher particle density and temperature. In pp collisions one 
expects no or very little collectivity; Fig. 9 indicates that the transition to this limit in the most peripheral 
Au+Au collisions is very steep, and that already a small number of nucleons participating in the collision 
generates significant collectivity and sizeable radial flow. 

Since the effects of the freeze-out temperature and flow on the spectral slope are strongly anticor- 
related [20], flow fits from single-particle spectra tend to produce narrow \ 2 valleys in the T-(v±) plane, 
oriented from the upper left to the lower right. A direct comparison of the spectra from Pb+Pb at the SPS 
and from Au+Au at RHIC shows that the RHIC spectra are flatter, especially for the heavier hadrons 
((anti)protons, (anti)Lambdas, etc.) [46, 57]. So, when fitted with exactly the same procedure, RHIC 
data give \ 2 valleys which in the T-(v±) plane lie to the upper right of those obtained from the SPS data 
[46, 57]. At the same Tf they produce a somewhat larger flow whereas with the same flow RHIC data 
yield somewhat higher freeze-out temperatures [46]. In either case the effect is about 15%. 

This flow analysis of the single-hadron spectra leads us to our First Lesson: When the collision 
is over and the hadron momentum spectra decouple, the fireball has a "thermal freeze-out temperature" 
Tf ~ 110 — 130 MeV and is in a state of rapid transverse collective expansion with (v±) ~ 0.45 — 0.55 c 
("The Little Bang"). These parameters have a systematic uncertainty ~ ± 10% and are strongly anti- 




correlated. In the T-(v±) plane the \ 2 valleys for the RHIC data are shifted systematically about 15% 
towards the upper right compared to the SPS data. 

4.1.9 Hydrodynamic calculation of radial flow 

Of course, a flow fit to the spectra is just a fit, and their is no a priori guarantee that the extracted fit pa- 
rameters make physical sense. In order to assess whether the extracted freeze-out temperatures and radial 
flow velocities are physically meaningful we must check that they can be the result of reasonable models 
for the dynamical evolution of the fireball. Is it possible to generate that much collective transverse flow 
during the short lifetime of the collision zone before everything fizzles? 

The natural language for studying flow phenomena is hydrodynamics. The hydrodynamic equa- 
tions of motion are obtained from the local conservation laws for energy, momentum and baryon number, 

d fl T^(x) = (i/ = 0,...,3) and d„j B {x) = 0, (32) 

by inserting the ideal fluid decompositions for the energy-momentum tensor and baryon number current: 

T^(x) = (e{x) +p(x)^u^{x)u u {x) - p(x) , f B {x) = n B {x)u^(x) . (33) 

These are 5 partial differential equations for 6 fields (p,e,riB and the 3 independent components of 
■u M ), and the system must be closed by providing an equation of state p(e, ub) for the fluid matter. 
We use an equation of state which mimicks the lattice QCD results from Fig. 2 by matching a hadron 
resonance gas below T cr to an ideal gas of massless quarks and gluons above T cr . For details I refer to the 
literature [58, 59]. The numerical solution is considerably simplified by assuming Bjorken scaling flow 
[41] along the beam direction [60] (see Sec. 4.1.4). The form of these simplified equations can be found 
e.g. in [58, 59]. They are solved for the transverse flow and hadron spectra near midrapidity where, at 
sufficiently high collision energies, the assumption of longitudinal boost-invariance should yield a good 
approxiumation. To obtain the entire rapidity distribution, on the other hand, requires a solution of the 
full (3+l)-dimensional set of equations (see, e.g., Refs. [61, 62]) since boost-invariance breaks down 
near the target and projectile rapidities. 

As discussed in Sec. 3.4, the applicability of fluid dynamics requires mean free paths which are 
short on any relevant macroscopic length scale. It can therefore not be used for the initial pre-equilibrium 
stage of the collision just after nuclear impact nor for the very last stage when the matter becomes so 
dilute that particles start to decouple. In addition to an equation of state (see Sec. 2.2), hydrodynamics 
thus needs initial conditions for the time at which it starts to be a good approximation and for the energy 
density, baryon number density and flow velocity profiles at that time, and & freeze-out algorithm which 
tells us where to stop the hydrodynamic evolution and how to translate the hydrodynamic output (energy 
and baryon densities and flow at freeze-out) into hadron spectra. For the freeze-out algorithm we use 
the Cooper-Frye formula discussed in Sec. 4.1.5 with a constant freeze-out temperature. The initial 
time, the initial entropy (or energy) and baryon density in the fireball, and the freeze-out temperature 
(or equivalently the decoupling energy density) are treated as parameters to be adjusted to experimental 
observables. The shape of the initial density profiles is calculated from the overlap geometry of the two 
nuclei (for details see [63]). For the initial flow we assume a Bjorken scaling profile along the beam 
direction (see Sec. 4.1.4) and zero initial transverse flow. 

We thus need to fix four parameters, r eq , s eq , n£ )Cq and ed ec - The ratio of the net baryon baryon 
density nB,eq to the entropy density s eq is fixed by the measured proton/pion ratio. By entropy conserva- 
tion, the final total charged multiplicity dN^/dy fixes the initial product (s • r) cq [50, 41, 60]. The value 
of r cq controls how much transverse flow can be generated until freeze-out. It can thus be determined 
by trying to fit the final pion and proton spectra simultaneously whose different slopes at low p± are 
sensitive to the radial flow at freeze-out. The associated freeze-out temperature, also implicit in these 
spectra, then fixes the decoupling energy density e^ec- 




Fig. 10: Identified pion, antiproton and kaon spectra for ^/snn = 130 GeV from the PHENIX [65, 69] and STAR [66, 67] 
collaborations in comparison with results from a hydrodynamic calculation [64]. The top left panel shows pion and (anti-)- 
proton spectra from central collisions. Shown in the other panels are spectra of five different centralities: from most central 
(top) to the most peripheral (bottom). The spectra are successively scaled by a factor 0.1 for clarity. 

The top left panel of Fig. 10 shows the hydrodynamic fit [64] to the transverse momentum spec- 
tra of positive pions and antiprotons, as measured by the PHENIX and STAR collaborations in central 
(b = 0) Au+Au collisions at -y/s = 130 A GeV [65, 66, 67]. The fit yields an initial central entropy den- 
sity s eq = 95 fm -3 at an equilibration time r eq = 0.6 fm. This corresponds to an initial temperature of 
T eq = 340 MeV and an initial energy density e = 25 GeV/fm 3 in the fireball center. (Note that these pa- 
rameters satisfy the "uncertainty relation" r eq • T eq « 1.) Freeze-out was implemented on a hypersurface 
of constant energy density with edec = 0.075 GeV/fm 3 . Table 1 summarizes the initial conditions applied 
in hydrodynamic studies at SPS and RHIC energies (see the recent review [58] for references). 
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17 
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43 
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T eq (MeV) 


257 


340 
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0.6 


0.6 



Table 1. Initial conditions for SPS and RHIC energies used to fit the particle spectra from central Pb+Pb or Au+Au collisions. 
Seq and T eq refer to the maximum values at r eq in the fireball center. 

The fit in the top left panel of Fig. 10 was performed with a chemical equilibrium equation of state. 
Use of such an equation of state implicitly assumes that even below the hadronization temperature T cr 
chemical equilibrium among the different hadron species can be maintained all the way down to kinetic 
freeze-out. With such an equation of state the decoupling energy ed ec = 0.075 GeV/fm 3 translates into a 
kinetic freeze-out temperature of Td ec ~ 130 MeV. The data, on the other hand, show [68] that the hadron 
abundances freeze out at T c h em ~ T cv , i.e. already when hadrons first coalesce from the expanding quark- 
gluon soup the inelastic processes which could transform different hadron species into each other are too 



slow to keep up with the expansion. The measured pj'K ratio thus does not agree with the one computed 
from the chemical equilibrium equation of state at the kinetic freeze-out temperature T& ec = 130 MeV, 
and the latter must be rescaled by hand if one wants to reproduce not only the shape, but also the correct 
normalization of the measured spectra in Fig. 10. 

A better procedure would be to use a chemical non-equilibrium equation of state for the hadronic 
phase [70, 7 1 , 72] in which for temperatures T below r c h em the chemical potentials for each hadronic 
species are readjusted in such a way that their total abundances (after decay of unstable resonances) are 
kept constant at the observed values. This approach has recently been applied [73] to newer RHIC data 
at ^ = 200 A GeV. 

Once the parameters have been fixed in central collisions, spectra at other centralities and for 
different hadron species can be predicted without introducing additional parameters. The remaining 
three panels of Fig. 10 show the transverse momentum spectra of pions, kaons and antiprotons in five 
different centrality bins as observed by the PHENIX [65, 69] and STAR [66, 67] collaborations. For 
all centrality classes, except the most peripheral one, the hydrodynamic predictions (solid lines) agree 
pretty well with the data. The kaon spectra are reproduced almost perfectly, but for pions the model 
consistently underpredicts the data at low p±. This has now been understood to be largely an artifact 
of having employed in these calculations a chemical equilibrium equation of state all the way down to 
kinetic freeze-out. More recent calculations [73] with a chemical non-equilibrium equation of state show 
that, as the system cools below the chemical freeze-out point T c h em ~ T cr , a significant positive pion 
chemical potential builds up, emphasizing the concave curvature of the spectrum from Bose effects and 
increasing the feeddown corrections from heavier resonances at low p±. 

Significant discrepancies are also seen at large impact parameters and large transverse momenta 
p± > 2.5GeV/c. This is not surprising since high-p^ particles require more rescatterings to thermalize 
and escape from the fireball before doing so. This is in particular true in more peripheral collisions where 
the reaction zone is smaller. 



Fig. 11: Transverse mass spectrum of f2 hyperons from central 
200 A GeV Au+Au collisions at RHIC [74]. The curves are hy- 
drodynamic calculations [73] with different initial and freeze- 
out conditions: Solid lines correspond to the default of no initial 
transverse flow at r cq , dashed lines assume a small but non-zero 
radial flow, v± = tanh(arx) with a = 0.02 fm _1 , already at 
T cq . The lower (thin) set of curves assumes f2-decoupling at 
T cr = 164 MeV, the upper (thick) set of curves decouples the £1 
together with the pions and protons at ed cc = 0.075 GeV/fm 3 . 
(In the chemical non-equilibrium hadronic equation of state em- 
ployed here this value of ed oc corresponds to a freeze-out tem- 
perature T dcc « 100 MeV - see Ref. [73] for details.) 

Figure 11 compares preliminary spectra of SI hyperons [74] with hydrodynamic predictions [73, 
75]. For this comparison the original calculations for 130 A GeV Au+Au collisions [75] were repeated 
with RHIC2 initial conditions and a chemical non-equilibrium equation of state in the hadronic phase 
[73]. Following a suggestion that U hyperons, being heavy and not having any known strong coupling 
resonances with pions, should not be able to participate in any increase of the radial flow during the 
hadronic phase and thus decouple early [76], Fig. 1 1 shows two solid lines, the steeper one corresponding 
to decoupling at e<j e c = 0.45 GeV/fm 3 , i.e. directly after hadronization at T cr , whereas the flatter one 
assumes decoupling together with pions and other hadrons at ed ec = 0.075 GeV/fm 3 . The data clearly 
favor the flatter curve, suggesting intense rescattering of the fi's in the hadronic phase. The microscopic 
mechanism for this rescattering is still unclear. However, without hadronic rescattering the hydrodynamic 




model, in spite of its perfect local thermalization during the early expansion stages, is unable to generate 
enough transverse flow to flatten the spectra as much as required by the data. Partonic hydrodynamic 
flow alone can not explain the Q spectrum. Similar comments apply to the Q and £1 spectra from Pb+Pb 
collisions at the SPS shown in Fig. 7. 

We close this subsection by concluding that the freeze-out parameters extracted from flow fits to 
the measured spectra are dynamically consistent, i.e. they can be reproduced in hydrodynamic calcula- 
tions of the fireball evolution using realistic initial conditions. The successful hydrodynamic calculations 
require fast thermalization at r eq < 1 fm/c and high initial energy density, e eq > 20 GeV/fm 3 . This con- 
clusion will be solidified by the study of elliptic flow in the next subsection. The hydrodynamically 
generated strong transverse flow naturally explains the much flatter proton than pion spectra at low p± 
(see Fig. 10). Together with the approximate baryon/antibaryon symmetry at RHIC (p/p « 0.7 in Au+Au 
collisions) this provides a simple understanding for the initially puzzling observation that in Au+Au col- 
lisions at RHIC for p± > 2 GeV/c antiprotons are as abundant as negative pions (upper left panel in 
Fig. 10)! This is completely unexpected from standard string and jet fragmentation phenomenology 
in pp and e + e~ collisions where the p/ir~ ratio never exceeds 10-15%, but a simple (albeit dramatic) 
consequence of the collective transverse flow in heavy-ion collisions. 



4.2 Elliptic flow 

4.2.1 Anisotropic flow in non-central collisions 

For central collisions (b = 0) between equal spherical nuclei, radial flow (left panel of Fig. 12) is the 
only possible type of transverse flow allowed by symmetry. In non-central (b / 0) collisions between 
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Fig. 12: Left and center panels: Schematic picture of radial and elliptic flow. Right panel: Initial energy density contours in the 
transverse plane for a Au+Au collision at impact parameter b = 7 fm. The dashed circles indicate transverse projections of the 
colliding nuclei. 

spherical nuclei (right panels of Fig. 12), or in b = collisions between suitably aligned deformed nuclei 
(such as 238 U), this azimuthal symmetry is broken and anisotropic transverse flow patterns can develop. 
The overlap region of the two colliding nuclei is then spatially deformed in the transverse plane, i.e. the 
spatial eccentricity of the reaction zone 
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(where the average is taken with the energy density distribution) is initially non-zero and positive. Rescat- 
tering processes among the produced particles transfer this spatial deformation onto momentum space, 
i.e. the initially locally isotropic transverse momentum distribution of the produced matter begins to 
become anisotropic. The driving force for this momentum anisotropy is the spatial eccentricity; the 



momentum anisotropy can only grow as long as e x > 0. As a function of time e x decreases, either spon- 
taneously due to free-streaming radial expansion (if no rescattering happens) or somewhat more quickly 
due to the development of elliptic flow (if rescattering occurs) which makes the system expand faster into 
the reaction plane than perpendicular to it (see middle panel in Fig. 12). Once the spatial eccentricity has 
disappeared, the momentum anisotropy saturates [77]. Since this happens quite early in the collision (as 
we will see shortly), the finally observed momentum anisotropy opens a window onto the early stage of 
the fireball expansion. 





Fig. 13: Snapshots of the energy density contours (lines) and velocity profiles (small arrows) taken at t— r oq = 3.2, 4.0, 5.6 and 
8.0 fm/c after initialization, for Pb+Pb collisions at b — 7 fm for an initial central temperature T cq = 500 MeV (from Ref. [59]). 
The outer thick (black) line indicates the freeze-out surface at edec, the next thick (blue) line separates the hadron resonance 
gas from the enclosed mixed phase, and the innermost thick (red) line separates the mixed phase from the quark-gluon plasma 
in the center. The quark-gluon plasma phase disappears shortly after the third snapshot at r— r eq = 5.6 fm/c (lower left panel). 
The horizontal and vertical axes are x and y in fm, respectively (as in Fig. 12). 

The mechanism for the creation of elliptic flow is most easily understood in the hydrodynamic 
limit: After thermalization there is high pressure in the interior of the reaction zone which falls off to 
zero outside. Obviously the pressure gradient is steeper in the short direction, leading to stronger hy- 
drodynamic acceleration into the reaction plane. This is seen in Fig. 13 where I show hydrodynamically 
computed energy density contours and flow vector fields for non-central Pb+Pb collisions for 4 different 
times after the onset of hydrodynamic expansion. One sees that the reaction zone expands faster into the 
reaction plane and that the spatial eccentricity is almost gone after about 8 fm/c. 

Figure 14 shows the time evolution of the spatial eccentricity e x and of the momentum anisotropy 

p{> fdxd y (T^ + Tyy)' 1 ; 

Note that with these sign conventions, the spatial eccentricity is positive for out-of-plane elongation (as 
is the case initially) whereas the momentum anisotropy is positive if the preferred flow direction is into 
the reaction plane. The calculations were done for Au+Au collisions at impact parameter 6 = 7 fm, for 



RHIC initial conditions with a realistic equation of state (EOS Q, solid lines) and for a much higher initial 
energy density (initial temperature at the fireball center = 2 GeV) with a massless ideal gas equation of 
state (EOS I, dashed lines), relevant for future experiments at the LHC [78]. 




Fig. 14: Time evolution of the spatial eccen- 
tricity e x and the momentum anisotropy e p for 
Au+Au collisions with b — 7 fm at RHIC and at 
much higher collision energies (from Ref. [78]). 



The initial spatial asymmetry at this impact parameter is e x (T cqu ) = 0.27, and clearly e p (T cqu ) = 
since the fluid is initially at rest in the transverse plane. The spatial eccentricity is seen to disappear 
before the fireball matter freezes out, in particular for the case with the very high initial temperature 
(dashed lines) where the source is seen to switch orientation after about 6 fm/c and becomes in-plane- 
elongated at late times [51]. One also sees that the momentum anisotropy e p saturates at about the same 
time when the spatial eccentricity e x vanishes; all of the momentum anisotropy is built up during the first 
5 fm/c. 

At RHIC the saturation happens even earlier because part of the matter begins to hadronize even 
before the spatial eccentricity has fully disappeared. Near the phase transition (in particular if as in our 
case it is modelled as a first order transition) the equation of state becomes very soft, and this inhibits 
the generation of transverse flow. This also affects the generation of transverse flow anisotropics as seen 
from the solid curves in Figure 14: The rapid initial rise of e v suddenly stops as a significant fraction 
of the fireball matter enters the mixed phase. It then even decreases somewhat as the system expands 
radially without further acceleration, thereby becoming more isotropic in both coordinate and momentum 
space. Only after the phase transition is complete and pressure gradients reappear, the system reacts to 
the remaining spatial eccentricity by a slight further increase of the momentum anisotropy. The softness 
of the equation of state near the phase transition thus focusses the generation of anisotropic flow to even 
earlier times, when the system is still entirely partonic and has not even begun to hadronize. At RHIC 
energies this means that almost all of the finally observed elliptic flow is created during the first 3-4 fm/c 
of the collision and reflects the hard QGP equation of state of an ideal gas of massless particles (c 2 s = |) 
[59]. Microscopic kinetic studies of the evolution of elliptic flow lead to similar estimates for this time 
scale [77,79, 80, 81]. 



4.2.2 The elliptic flow coefficient V2 

The momentum anisotropy e p manifests itself as an azimuthal anisotropy of the measured hadron spectra. 
One quantifies this anisotropy in terms of the azimuthal Fourier coefficients of the transverse momentum 
spectrum: 

These coefficients depend on the impact parameter b, the transverse momentum p± and the particle 
species i (through their rest mass mi). We have suppressed the dependence on rapidity y since the cal- 
culation assumes longitudinal boost-invariance and thus makes meaningful statements only near y = 0. 
At midrapidity V2 is (for collisions between equal nuclei) the lowest nonvanishing Fourier coefficient. 



Reflection symmetry with respect to the reaction plane forbids the appearance of sine terms in the ex- 
pansion. Equation (36) implies that the elliptic flow coefficient can be calculated as 

4(p±,6) = (cob(2^))* Xj6 (37) 

were the average is performed with the transverse momentum spectrum of particle species i at fixed 
impact parameter and p±. If the average is performed with the p± -integrated spectrum one obtains the 
"P± -integrated elliptic flow" v l 2 (b). 




Fig. 15: Elliptic flow of all charged hadrons in 130/lGeV Au+Au collisions at RHIC. Left: p± -integrated elliptic flow vs. 
collision centrality (measured by the normalized total charged particle multiplicity n c h). Right: differential elliptic flow V2(p±) 
for minimum bias collisions (i.e. inegrated over all impact parameters). The curves are hydrodynamic calculations for various 
parametrizations of the initial transverse density profile (for details see [63]). Data are from the STAR [82] and PHENIX [83] 
Collaborations. 

Figure 15 shows a comparison of elliptic flow data for charged hadrons from Au+Au collisions 
at RHIC with hydrodynamic calculations. The hydrodynamic results are genuine predictions since all 
model parameters have been fixed in central (b = 0) collisions where V2 vanishes. The impact param- 
eter dependence of V2(p±) then follows from the overlap geometry and its implications for the initial 
transverse density profile and its spatial eccentricity [63], without additional parameters. The data follow 
the hydrodynamic predictions up to impact parameters b > 7fm (n c h/n max < 0.5 in Fig. 15a) and up to 
transverse momenta pj_ < 2GeV/c (Fig. 15b). For p± > 2GeV/c (note that this concerns less than 0.5% 
of all produced hadrons!) the measured elliptic flow saturates and falls behind the hydrodynamically 
predicted continuous rise with p±. This indicates a gradual breakdown of local thermal equilibrium for 
high-p^ particles - not a big surprise. Similarly, the elliptic flow in peripheral collisions is also smaller 
than predicted - again presumably a consequence of incomplete local thermalization in very peripheral 
collisions where the nuclear overlap zone becomes quite small. 

Hydrodynamics predicts a clear mass-ordering of elliptic flow [75]. As the collective radial motion 
boosts particles to higher average velocities, heavier particles gain more momentum than lighter ones, 
leading to the previously discussed flattening of their spectra at low transverse kinetic energies. When 
plotted against p± this effect is further enhanced by a kinematic factor arising from the transformation 
from m± to p± (see Fig. 10). This flattening reduces the momentum anisotropy coefficient v-i at low p± 
[75], and the heavier the particle the more the rise of V2(p±) is shifted towards larger p±. 

This effect, which is a consequence of both the thermal shape of the single-particle spectra at low 
p± and the superimposed collective radial flow, has been nicely confirmed by the experiments: Figure 16 
shows that the data [84, 85] follow the predicted mass ordering out to transverse momenta of about 
1.5 GeV/c. For K® and A+A very accurate data have recently become available from 200 AGeV Au+Au 
collisions [86], again in quantitative agreement with hydrodynamic calculations up to p± ~ 1.5 GeV/c for 




Fig. 16: Differential elliptic flow for identified pions, kaons and (anti)protons from Au+Au collisions at ^/snn = 130 GeV [84] 
(left) and at ^/snn = 130 GeV [85] (right), together with hydrodynamic calculations [64]. 

kaons and up to p± ~ 2.5GeV/c for A + A. The inversion of the mass-ordering in the data at large pj_ 
is caused by the mesons whose V2(p±) breaks away from the hydrodynamic rise and begins to saturate 
at p± > 1.5 GeV/c. In contrast, baryons appear to behave hydrodynamically to p± ~ 2.5 GeV/c, breaking 
away from the flow prediction and saturating at significantly larger pj_ than the mesons. This is consis- 
tent with the idea that the hadronic elleptic flow really reflects a partonic elliptic flow already established 
before hadronization, that the latter exhibits a hydrodynamic rise at low pj_ followed by saturation above 
p± ~ 750 — 800 MeV/c, and that these features are transferred to the observed hadrons by quark coales- 
cence, manifesting themselves there at twice resp. three times larger p± -values [87]. 

The comparison of the V2 data with the hydrodynamic model makes the implicit assumption that all 
of the measured V2 is collective, i.e. that non-flow contributions to the angular average V2 = (cos(2<p p )) 
can be eliminated. Various experimental cuts [82] ensure that the influence of non-collective 2-particle 
correlation effects (jet correlations [89], resonance decays, HBT effects) on V2 are minimized. An alter- 
native method to extract V2 from higher-order cumulants of the azimuthal particle distribution [90, 91] 
has recently shown [92] that, in the interesting range p± < 2 GeV/c where the emitted hadrons behave 
hydrodynamically, non-flow corrections do not exceed 15% of the value extracted from the prescription 
V2 = (cos(2ip p )}. (This may be different at much higher pj_ where V2 from 2-particle angular correla- 
tions seems to approach a p± -independent constant value [93] but statistics is not good enough for a 
4-particle cumulant analysis which could exclude that at high p± an increasing fraction of this V2 stems 
from angular 2-particle jet correlations.) Further support of the collective nature of the measured V2 
comes from the observation that the event plane angles reconstructed from V2 at forward and backward 
rapidities are, within the statistical error, perfectly correlated. We can therefore assume that at least 85% 
of the measured V2 values shown in Figs. 15-17 are due to collective flow. 

4.2.3 Rapid thermalization and early pressure — the creation of the Quark-Gluon Plasma 

I already mentioned that momentum anisotropics, in particular elliptic flow, can only build as long as the 
source is spatially deformed. Even if the created particles don't rescatter but simply fly off ballistically 
with their isotropically distributed initial momenta, the spatial eccentricity decreases as a function of 
time as [59] 

s x (t +At) 

where to is the time when the particles were created and (f 2 ) T0 is the azimuthally averaged initial 
transverse radius squared of the reaction zone. So if thermalization is delayed by a time At, any elliptic 
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flow would have to build on a reduced spatial deformation and would come out smaller. 

It is also known from microscopic kinetic studies [80, 81] that, for a given collision geometry, i.e. 
a given initial spatial eccentricity e x , the magnitude of the generated elliptic flow is a monotonic function 
of the mean free path (or the product of density and scattering cross section) in the fireball. This is shown 
in Fig. 17. Without reinteractions (i.e. for infinite mean free path) no momentum anisotropics develop at 
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all [80]; the maximal value for V2 is reached in the hydrodynamic limit of zero mean free path [88]. 

This hydrodynamic upper limit for the charged particle elliptic flow was shown to be insensitive to 
shape details of the initial transverse density distribution [63] and, as long as the shape and normalization 
of the single particle spectra from central collisions were held fixed by proper retuning of the initial 
conditions, any sensitivity to details of the QGP equation of state remained below 25% [75]. (Note that 
this statement refers to the V2 of all charged particles which are strongly dominated by pions and kaons; 
Fig. 16a shows that the elliptic flow of heavier particles such as protons, which are more sensitive to 
radial flow, shows indeed significant dependence on the equation of state, and that the data disfavor, for 
example, a pure hadron resonance gas over an equation of state with a deconfining phase transition and 
a hard QGP equation of state above T cr .) 

The observation that for b < 8 fm and pj_ < 1.5 — 2 GeV/c the measured V2 exhausts at least 85% of 
the hydrodynamically predicted value thus puts severe constraints on the dynamical evolution and on our 
global picture of the collision event. It requires early thermalization (such that the spatial eccentricity 
does not decay before elliptic flow starts to build) and almost complete local equilibrium throughout the 
crucial first 5 fm/c of the expansion when the momentum anisotropics are generated. The hydrodynamic 
model calculations for RHIC were done with a thermalization time r eq = 0.6 fm/c; increasing r cq beyond 
1 fm/c begins to not only seriously degrade the description of the angle-averaged single-particle spectra, 
but also cuts down v-i below the experimentally measured level. To obtain the full measured amount of 
elliptic flow requires almost complete thermalization in less than 1 fm/c. 

The sensitivity of V2 to any kind of deviation from local thermal equilibrium seems to be particu- 
larly strong [94, 95]. (Of course, we are talking here about the momentum distribution for particles with 
average momenta p± ~ 3 T; no claim for thermalization can be made for particles with much larger than 
average p± because they are too rare to affect bulk properties such as the pressure and the hydrodynamic 
expansion.) Turning this observation around and remembering that most of the momentum anisotropy is 
generated at energy densities well above the critical value for hadronization, the good agreement of the 
data with ideal fluid dynamics points to a very small viscosity of the quark-gluon plasma. (Note my use 
of the word quark-gluon plasma for the matter in the early expansion stage! What other name could I use 
for a thermalized system of strongly interacting matter at energy densities > 10 e cr ?) The QGP does not 
behave as a weakly interacting quark-gluon gas, as suggested by naive perturbation theory, nor does it 



behave like viscous honey, as suggested by the value of the shear viscosity extracted from state-of-the-art 
resummed thermal field theory at leading order in the strong coupling constant [36]. The fact that the 
QGP behaves like an ideal fluid implies strong non-perturbative interactions in the quark-gluon plasma 
phase. 

So here is our Second Lesson: The strong measured elliptic flow at RHIC requires rapid thermal- 
ization at r cq < 1 fm/c. At this time the average energy density in the thermalized fireball is e > 10— 15e cr , 
so thermalization leads into a quark-gluon plasma state. In Au+Au collisions at RHIC this QGP lives for 
about 5-7 fm/c before hadronizing. 

So what about the SPS? Have we made quark-gluon plasma there, too [10]? The energy densities 
e > 3.5GeV/fm 3 were sufficient and, on the 30% accuracy level, hydrodynamics did a good job for the 
angle-averaged single particle spectra from S+S and Pb+Pb collisions at the SPS, too [20, 96]. However, 
the measured elliptic flow at the SPS is only about 50% of the value measured at RHIC [97, 98] whereas 
hydrodynamics predicts it to be even slightly larger than at RHIC [59]. 1 So thermalization seems to be 
less efficient at the SPS than at RHIC. This may not be unreasonable, given the significantly lower (by 
about a factor 3 [59]) initial density and the rapid conversion of the matter to a less strongly coupled 
hadron resonance gas before the momentum anisotropy has been able to reach saturation. It would be 
good to test this with central U+U collisions at the lower end of the collision energy range of RHIC 
[59]. Such collisions provide similar spatial eccentricities e x as semiperipheral Au+Au collisions, but 
for a fireball which covers a much larger transverse area. This should improve thermalization also at 
lower energies, making hydrodynamics a better approximation, and help to test our understanding of 
the dynamical origin of elliptic flow. - I should also mention that at RHIC hydrodynamics fails to 
reproduce the rapidity-dependence of elliptic flow, predicting a wide rapidity plateau for V2 [72] whereas 
the measured elliptic flow rapidly decreases as one moves away from midrapidity [99]. The origin of 
this discrepancy is not yet clear, but I suspect that it again indicates a breakdown of local equilibrium 
away from midrapidity [100], for reasons which are related to the analogous breakdown at midrapidity 
at lower collision energies [101]. 

5 TWO-PARTICLE BOSE-EINSTEIN CORRELATIONS AND HBT INTERFEROMETRY 

Single-particle momentum spectra, including their anisotropics, provide no direct information on the 
space-time structure of the source. It is therefore in principle possible to describe these spectra with 
models that have the same momentum structure but completely different freeze-out distributions in space 
and time. Two-particle correlations, on the other hand, depend on the average separation of the particles 
at decoupling and therefore provide valuable spatial and temporal information. This is most easily seen 
for correlations generated by final state Coulomb interactions between charged hadrons which influence 
their distribution of relative momenta depending on how close they were at time of emission. The 
measured relative momentum distribution of particle pairs thus allows to constrain their average space- 
time distance at decoupling. Conceptually more involved, but theoretically easier to analyze are two- 
particle momentum correlations between pairs of identical particles ("Bose-Einstein correlations" in the 
case of identical bosons) caused by quantum statistical (wave function symmetrization) effects. Again, 
these affect the relative momentum distribution at low relative momenta, and the relative momentum 
range of these quantum statistical correlations yields information about the average space-time separation 
of the emission points of the two particles, and thereby also about the overall size and emission duration 
of the source. There exist a number recent reviews of these methods to which I refer the reader for more 
details [102, 103, 104, 105]. 

For pairs of identical bosons (e.g. pions), this technique is known as Hanbury Brown-Twiss (HBT) 
interferometry [106]. The basic principle, as applied in particle and nuclear physics, is illustrated in 

'At the SPS most of the elliptic flow is generated in the hadronic phase, and the softening effects of the phase transition on 
the development of flow are not as important as at RHIC [59]. 




Fig. 18: The principle of two-particle interferometry: a pair of indistinguishable particles is emitted from points a and b and 
detected by detectors 1 and 2 which measure their momenta pi and p2- The total pair yield is obtained by integrating over all 
possible positions (a, b) in the source. The measurement is performed as a function of the average momentum K — (p 1 +P2 ) /2 
of the pair and its relative momentum q = pi — p2- The q-dependence yields information about the average separation (R 2 ) 1 ^ 2 
of the emitters a and b which contribute pairs of given pair momentum K . 

Figure 18. I will now shortly outline its theoretical treatment, referring to the literature [102, 103, 104, 
105] for more details. 



5.1 Single-particle spectra and two-particle correlations 

The covariant single- and two-particle distributions are defined by 

Pi(p) = E^- = E(a+a p ), (39) 

Pi (Pi i P2 ) = El E2 d 3p^ 3p2 = E l E 2 («p! «p 2 «P2 Opi ) , (40) 

where a+ (d p ) creates (destroys) a particle with momentum p. The angular brackets denote an ensemble 
average (O) = tr (pO) where p is the density operator associated with the ensemble. (When talking 
about an ensemble we may think of either a single large, thermalized source, or a large number of 
similar, but not necessarily thermalized collision events.) The single-particle spectrum is normalized to 
the average number of particles, (N), per collision, 

J^P 1 (p) = (N), (41) 

while the two-particle distribution is normalized to the number of particles in pairs, (N(N — 1)), per 
event: 

' /> d> '' 2 P2(px,P2) = (N(N-l)). (42) 
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The two-particle correlation function is defined as 

C{Pl,P2) = p \ r) / \ ■ ( 43 ) 

If the two particles are emitted independently and final state interactions are neglected one can prove a 
generalized Wick theorem [107, 108] 

|(ajr a„ )| 2 

C(pi,P2) = 1 ± : ■ (44) 

{a Pl a Pl ){ap 2 a P2 ) 



This assumption is often called "chaoticity assumption", having in mind a completely chaotic (uncorre- 
lated) emitter, in contrast to processes where the phases of the wave functions of the two emitted particles 
are to some extent correlated with each other. The degree of "chaoticity" or "phase correlation" in the 
source can be studied by comparing three-particle correlations with two-particle ones [109]. This was 
recently done by the STAR Collaboration [110] who established that in central Au+Au collisions at 
a/ s nn = 130 GeV the pion emitting source is indeed completely chaotic. 

From now on I will assume that the emitted particles are bosons, e.g. pions, and use the + sign in 
Eq. (44). We next discuss the numerator of the Bose-Einstein correlation term, Ka^cip^] 2 . 



5.2 Source Wigner function and spectra 

As reviewed in Refs. [103, 107, 108], the pion emitting source can be described in terms of classical 
currents J(x) which act as classical sources of freely propagating pions. These currents represent a 
parametrization of the last collision from which the free outgoing pion emerges. Using this language in 
one form or another [103, 107, 108], one can express the "exchange amplitude" (a^a^) as the Fourier 
transform of the single-particle Wigner density S(x, K) which is the quantum mechanical analogue of 
the classical phase-space density at pion freeze-out, describing the probability of emitting a pion with 
momentum K from space-time point x: 



J d 4 xS(x,K) e^ 1 



(45) 



S(x, K) is also called "emission function". K = ^(Ei+E2,pi+P2) is the average momentum of the 
pair while q = (Ei—E 2 ,Pi—p2) is the relative momentum between the two pions. Since the pion 
momenta p\ , p 2 are on-shell, = Ei = (m 2 + p 2 ) 1 ^ 2 , the 4-momenta q and K are in general off-shell. 
They satisfy the orthogonality relation ("on-shell constraint") 



q-K = 



q° = (3 ■ q with 
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(46) 



(3 is the velocity of the pair; the last, so-called "on-shell" approximation holds for small relative momenta 
q which is where the correlation term is non-zero: 




(47) 



Since the relevant range of q is given by the inverse size of the source (more properly: the inverse size 
of the regions of homogeneity in the source - see below), the validity of this approximation is ensured in 
practice as long as the Compton wavelength of the particles is small compared to this "source size". For 
the case of pion, kaon, or proton interferometry for heavy-ion collisions this is true automatically due 
to the rest mass of the particles: even for pions at rest, the Compton wavelength of 1.4 fm is comfort- 
ably smaller than any typical nuclear source size. This is of enormous practical importance because it 
allows you essentially to replace the source Wigner density by a classical phase-space distribution func- 
tion for on-shell particles. This provides a necessary theoretical foundation for the calculation of HBT 
correlations from classical hydrodynamic or kinetic (e.g. cascade) simulations of the collision. 

With the definition (45) we can express the single-particle spectra (39) and correlation function 
(44) through the emission function S(x, K): 

^ d A x S(x,p) , 
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The last line uses the "smoothness approximation" P\(K+^q)Pi(K— |qr) « \P\(K)\ 2 which can be 
justified for sufficiently large sources; it is not really necessary because the single-particle spectra are 
measured and the ratio \P\(K)\ 2 / Pi(K+^q)Pi(K—^q) can thus be divided out from the correlator 
C(q, K) — 1 before analyzing it theoretically, but we use it here for simplicity. 

The last equality in Eq. (49) introduces the notation (f(x)) for the average of a space-time ob- 
servable f(x) with the emission function S(x, K). Note that, due to the independence of the emission 
function, this average depends on the pair momentum K. 

The fundamental relations (48) and (49) show that both the single-particle spectrum and the two- 
particle correlation function can be expressed as simple integrals over the emission function. Whereas 
the single-particle momentum spectrum (48) integrates over the spatial structure of the emission function 
and thus provides no spatial information about the source, Eq. (49) shows that the g-dependence of the 
correlation function C(q, K) provides access to the space-time structure of the source. However, it de- 
pends only on the modulus square of the spatial Fourier transform of the source, so the phase information 
is missing, and furthermore, due to the "on-shell constraint" (46), not all four components of q can be 
varied independently. The Fourier transform in Eq. (49) thus mixes space and time in a particular way, 
and the emission function S(x, K) cannot be fully reconstructed from the measured correlation func- 
tion C(q, K) without taking recourse to model assumptions about the source. Still, as I will now show, 
the correlator C(q, K) provides very important constraints for any dynamical model for the heavy-ion 
collision. 

5.3 HBT radii and source sizes 

The most interesting feature of the two-particle correlation function is its half-width as a function of 
q. Actually, since the relative momentum q = p\—p2 has three Cartesian components, the fall-off of 
the correlator for increasing q is not described by a single half- width, but rather by a (symmetric) 3x3 
tensor. Through the Fourier transform in Eq. (49), these half-widths of the correlator (at a given pair 
momentum K) can be related to the widths of the emission function S(x, K) in space and time (at the 
same K). This relation becomes exact if particle emission is not complicated by processes involving 
multiple time scales. Multiple time scales, with corresponding dynamically induced multiple length 
scales, are an issue in pion emission from heavy-ion collisions: a large fraction of the observed pions is 
not emitted from the source directly, but appears later from the decay of unstable resonances which were 
originally emitted directly from the source but than travelled a certain distance without interactions before 
decaying. Resonance decay effects on the HBT correlation function have been extensively discussed in 
the literature (see e.g. [Ill, 112]). The results can be shortly summarized: Most resonances are so 
short-lived that they don't travel far before decaying, leading to a pion emission function of almost the 
same size and shape as for directly emitted pions. Some resonances are very long-lived and decay far 
outside the original source; the corresponding temporal and spatial tails in the emission function affect 
the correlator at such small values of q that they can not be experimentally resolved. They effectively 
reduce the visible Bose-Einstein correlation strength and can be absorbed by reducing the second term in 
Eq. (49) by a K-dependent "intercept parameter" \{K): C(q,K) = 1 + X(K)\ {e iq ' x ) dil \ 2 (K) where 
(. . .)dir denotes the average with the emission function for (almost) directly emitted pions. The only 
resonance that can neither be classified as short-lived nor as long-lived in this sense is the oj meson with 
its lifetime of about 20fm/c. I will assume here that its abundance is small enough to be negligible 
(this depends on the freeze-out temperature), but this is mostly for didactical purposes and needs to be 
reassessed in the actual comparison between theory and data. 

With these assumptions the "direct" emission function can be fully characterized by its (in- 
dependent) "widths" 

Sn„(K) = [(x^)^ - (x M ) dir (x i ,) dir )] (K) = (x^x v ) dir (50) 
which describe the widths of the effective emission region inside the source ("homogeneity region") 



which effectively contributes to the emission of particles with momentum K . Note that for an expanding 
source these effective emission regions do not coincide with the entire source: As illustrated in Fig. 19, 
pions whose momentum K is large and points in the positive x direction will predominantly originate 
from the elliptic region near the right edge of the figure, since this region anyway moves collectively 
in the +x direction. On the other hand, the elliptic region near the top of the figure moves in the +y 
direction and thus mostly contributes pions with (large) positive y-momentum K y . Eq. (49) shows that 
the correlation function C(q, K) of pairs with pair momentum K probes only these effective emission 
regions for particles with fixed K, not the entire source. Only for static fireball, where each volume 
element radiates the same momentum spectrum, does the HBT correlation function measure the entire 
source. (The notation x s ;d e and x out used in Fig. 18 will be explained below.) 



If, as assumed, the emission function S(x, K) can be fully characterized by its "spatial correlation 
tensor" S^ U (K) (50) which describes the size and shape of the effective emission regions for particles 
of momentum K, we can make the following Gaussian ansatz for S(x, K): 



The normalization is fixed by the single-particle spectrum [103] but drops out from the correlation func- 
tion. Note the dependence only on the distance x^{K) = x^ — {x fl ){K) to the point of maximum 
emissivity at momentum K - the latter itself can not be measured, i.e. the absolute position of the 
regions of homogeneity is experimentally inaccessible. Reconstructing the entire source as a sum of 
effective emission regions for fixed K is therefore always a model-dependent theoretical enterprise. 

The Fourier transform of the above Gaussian parametrization is easily evaluated. Eliminating 
q° with the on-shell constraint (46) and introducing the intercept parameter X(K) as described above 
to account for the unresolvable tail in the emission function from long-lived resonances gives for the 
correlation function the Gaussian form 




Fig. 19: Illustration of the concept of homogeneity regions in an expanding source (see text). 
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(53) 



The Rij(K) are called "HBT radius parameters"; note that the off-diagonal terms with i ^ j need not be 
positive. The last equation shows explicitly the mixing between spatial and temporal source information 
in the experimentally measured HBT radius parameters, caused by the on-shell constraint (46). 

Note that this mixing involves the components of the pair velocity (3; it can be minimized by 
rotating the Cartesian coordinate system (x,y,z) such that K and (3 point in +x' direction [113]. In 
relativistic heavy-ion collisions it is, however, advantageous to keep the z-axis aligned with the beam di- 
rection, and therefore one rotates only the tranverse (x, y) axes such that the transverse pair momentum 
K± points in x' direction. This eliminates the y' components of K± and (3. Traditionally [113], the ro- 
tated x' axis is called the "outward" direction (denoted as x out or x Q ), the y' axis is called the "sideward" 
direction (x si( jc or x s ), while the beam direction z is called "longitudinal" (denoted by xi ong or xi). In 
this frame [3 = ((3±,0,(3l) and K = (K±,0, Kl)- Note, however, that in non-central heavy-ion col- 
lisions the direction of the impact parameter b defines a distinguished direction in the transverse plane, 
and that therefore even in the out-side-long system the HBT radius parameters Rij(K) still depend on 
both magnitude and direction of K±. 



5.4 Central (6=0) heavy-ion collisions 

5.4.1 Formalism 

In central collisions between spherical nuclei, the entire source is azimuthally symmetric around the 
beam axis. After fixing the transverse coordinate system such that the "out"-axis points in the emission 
direction K± of the pair, the HBT radii Rij can only depend on the magnitude, but not the direction of 
K±. Furthermore, the source is clearly symmetric under reflections with respect to the plane defined by 
the beam direction and K±, so all components of the spatial correlation tensor S^ v = {x^x u ) which are 
linear in x s id e must vanish. As a result, only six nonvanishing components of S^(K±, Y) survive. (I 
here introduce the pair rapidity Y = tanh^ 1 0l instead of the longitudinal pair momentum K^.) On the 
other hand, the same reflection symmetry also eliminates terms from the exponent in Eq. (52) that are 
linear in q s = q s idc such that only four HBT radius parameters are experimentally accessible: 



C(q, K) = 1 + \(K) exp -q 2 s R 2 s - q 2 R 2 - q 2 R 2 - 2q qiR 2 ol 



(54) 



The HBT radii are again functions of two variables only, R 2 j(K±, Y). Still, there are two more functions 
Sfj,v(K±, Y) needed to fully characterize the source than can be measured, due to the on-shell constraint. 

In symmetric A + A collisions at midrapidity Y = 0, things simplify further: the cross-term R 2 ol 
vanishes [114], leaving only three functions R S (K±), R (K±), and Ri(K±) to be determined. They are 
related to the source widths = {x^x u ) as follows: 

R 2 = (x 2 ), (55) 
R 2 Q = (x 2 )-2(3 ± (x t) + (3 2 ± (P), (56) 
Rf = (z 2 ), (57) 

where I used that 0i = at Y = 0. I will concentrate on this simple limit since RHIC HBT data have so 
far only been obtained near Y = 0. 

The HBT radii (55-57) depend on K± both explicitly (through the (3± factors in (56)) and implicitly 
(since (...) denotes an average with the iv"-dependent emission function). If K± is the only vector in 
the problem, there is no distinction between the outward and sideward directions in the limit K± — > 0, so 
R Q and R s must be equal to each other in that limit: 

i R o- R ") = K l % " = °- (58) 

This theorem is not valid for non-central collisions, where the impact parameter b introduces a preferred 
direction, and it may also break down for "opaque" (surface emitting) sources where, at least in principle, 
the outward normal vector can play a similar role [1 15, 1 16]. 



The geometric contributions (ot 2 ), {x 2 ), and (z 2 ) are controlled by an interplay between thermal 
motion and collective expansion: stronger expansion, resulting in larger flow velocity gradients, tends to 
reduce these variances whereas higher freeze-out temperatures, resulting in broader momentum spectra, 
tend to increase the size of the homogeneity regions. This size increase due to thermal smearing is easily 
understood by considering the limit T — > 0: in this limit particles have no thermal motion, and particles 
of momentum K can therefore only come from individual points in the source whose collective flow 
velocity agrees exactly with the pair velocity (3. The corresponding homogeneity regions are pointlike, 
i.e. have zero size. Thermal particle motion can smear out the collective flow velocity gradients such 
that particles of given K can now come from a larger homogeneity region. As T — > oo, the homegeneity 
regions approach the size of the entire fireball, irrespective of its dynamical state. 

Making these features explicit requires a parametrization of the emission function, for example 
the extensively studied model function [102, 103, 114, 117] 
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(59) 



Here the transverse radius r 2 = x 2 +y 2 , the spacetime rapidity i] = \ \n[(t+z)/(t— z)\, and the lon- 
gitudinal proper time r = Vt 2 —z 2 parametrize the spacetime coordinates with measure d A x = 
rdrdrirdrdip. Y = \ \n[(E K +K L ) / {E K -K L )\ and M ± = ^m 2 +K\ parametrize the longitudinal 

and transverse components of the pair momentum K. \[2R is the transverse geometric (Gaussian) ra- 
dius of the source, tq its average freeze-out proper time, At the mean proper time duration of particle 
emission, and A77 parametrizes the finite longitudinal extension of the source. T(x) is the freeze-out 
temperature; if you don't like the idea of invoking thermal equilibrium during the decoupling stage, you 
can think of T as a parameter that describes the random distribution of the particle momenta at each 
space-time point around their average value. The latter is parametrized by a collective flow velocity 
u^{x) in the form 

u^(x) = (cosh 77 cosh r]t(r), sinh r\ t (j) e r , sinh 77 cosh rjt(r)) , (60) 

with a boost-invariant longitudinal flow rapidity 77; = 77 (vi = z/t) and a linear transverse flow rapidity 
profile 

Vt(r) = Vf (j^j ■ (61) 

r/f scales the strength of the transverse flow. The exponent of the Boltzmann factor in (59) can then be 
written as 

K ■ u(x) = Mj_ cosh(Y — 77) cosh77t(r) — K±-e r sinh^r) . (62) 

For vanishing transverse flow (77/ = 0) the source depends only on M± and remains azimuthally sym- 
metric for all K±. In this case the sideward and longitudinal HBT radii (55) and (57) exhibit perfect 
Mi-scaling, i.e. when plotted as functions of Mj_ they coincide for pion and kaon pairs [103]. [In fact, 
R s does not depend on M± at all in this limit.] For the outward radius Mi-scaling is broken by the 
/3i -dependent terms in Eq. (56). For non-zero transverse flow, R s does depend on M±, at a level that 
directly reflects the strength of the transverse flow (see below). In this case the M^-scaling of R s and Ri 
is broken by the K± -dependent second term in the thermal exponent (62), albeit only weakly [103]. 

Using saddle-point integration, the HBT radii (55-57) for this source can be approximately evalu- 
ated analytically [1 14]. For pairs with Y = one finds the simple pocket formulae 

R 2 = Rl, (63) 
R 2 o = Rl + (3 2 ± (AQ 2 , (64) 
R? = Ll, (65) 



with 



1 1 1 

IV = i?2 + R^' (66) 

n * n -"flow 



(At,) 2 = (Ar) 2 + 2( v /r 2 + L2-r ) 2 , 



(67) 



1 11 

n = J^V 2 + ^' (68) 

Here Ra ow and Lfl ow are the transverse and longitudinal "dynamical lengths of homogeneity" generated 
by the expansion velocity gradients: 

£W(M ± ) = ^/^ = 75 _ 7 A_ t /__, (69) 




£flow(M ± ) = To J— = —J — , (70) 



where is the longitudinal 4-velocity. These expressions show explicitly the competition between flow 
gradients and thermal smearing (reflected in the factor y/T/M±) which I mentioned above. They also 
show the competition between the overall geometric size of the source and the "dynamical homogeneity 
lengths", with the smaller of the two controlling the HBT radii. The second term in the expression (67) for 
the emission duration reflects the fact that the emission duration is measured in terms of the coordinate 
time in the fixed center-of-mass reference frame, whereas freeze-out occurs at constant proper time tq 
which, over the longitudinal range probed by the longitudinal HBT radius, explores a finite range of 
coordinate times even if the proper time spread At is set to zero. For freeze-out along a hyperbola of 
fixed proper time, it is therefore impossible to get a vanishing emission duration At*. 

For a longitudinally boost-invariant source (A77 — ► 00) and weak transverse expansion the source 
function integral in Eq. (57) can actually be done exactly [118], and one finds instead of (65,70) 

02 r 2 T K 2 (MJT) 
For small M± ~ T, this differs by up to a factor of 2 from Eq. (70). 



5.4.2 HBT radii from Au+Au collisions at RHIC — the "HBT Puzzle " 

Let us now study the HBT radii predicted by the hydrodynamic model introduced in Section 4.1.9, which 
was very successful in reproducing the single particle momentum spectra, and compare them with data 
from RHIC. 

The hydrodynamic emission function, which implements sudden decoupling on a sharp freeze-out 
surface £ as described in Section 4.1.5, takes the form [119, 120] 

gi f K-d 3 a{x')S 4 (x-x'^ 



K) (2vr)3 h exp{[K-u(x>) - ^(x')}/T dec (x')} ± 1 ' (72) 

Phenomenological fits to spectra and HBT data often use a generalization of this form which replaces 
the ^-function by allowing for a spread of emission times ("fuzzy freeze-out") [114, 117, 121]. 

For a longitudinally boost invariant source the freeze-out hypersurface can be parametrized in 
terms of the freeze-out eigentime as a function of the transverse coordinates, Tf(x,y). The normal 
vector <i 3 <7 M on such a surface is given by 

d 3 a = ^cosfiT?, Vj_T/(a;,y), sinhr]jTf(x,y)dxdydr]. (73) 



With the four momentum = (M± cosh Y, K±, M± smh Y), where Y and M± are the rapidity and 
transverse mass associated with K, Eq. (72) becomes 



Si(x,K) 
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(2tt) s 



drj dx dy[M± cosh(y — rj)—K± • Vj_r/(;r, y)] 
xf(K-u(x),x) 5(r-Tf(x,y)) 



(74) 



with the flow-boosted local equilibrium distribution f(K-u(x),x) from (10). 

With this expression we can now study the emissivity of the source as a function of mass and 
momentum of the particles. For the purpose of presentation we integrate the emission function over two 
of the four space-time coordinates and discuss the contours of equal emission density in the remaining 
two coordinates. We begin with calculations describing central Au+Au collisions at ^/snn = 130 GeV 
and focus on directly emitted pions, neglecting pions from unstable resonance decays. Resonance decay 
pions are known to produce non-Gaussian tails in the spatial emission distribution, increasing its width, 
but these tails are not efficiently picked up by a Gaussian fit to the width of the measured two-particle 
momentum correlation function [112, 122]. A comparison of the experimental HBT size parameters 
extracted from such fits with the spatial widths of the emission function is thus best performed by plotting 
the latter without resonance decay contributions. 
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Fig. 20: Pion source function S(x, K) for central Au+Au collisions at ^/ snn = 130 GeV. The upper row shows the source after 
integrating out the longitudinal and temporal coordinate, in the lower row the source is integrated over the longitudinal and one 
transverse coordinate (y). In the left column we investigate the case K = 0, in the right column the pions have rapidity Y = 
and transverse momentum K± = 0.5 GeV in x direction [88]. 

Figure 20 shows equal density contours at 10, 20, . . . , 90% of the maximum in a transverse cut 
integrated over time and rj (top row) and as a function of radius and time integrated over r] and the 
second transverse coordinate (bottom row). The dashed circle in the top row indicates the largest freeze- 
out radius reached during the expansion, the dashed line in the bottom row gives the freeze-out surface 
Tf(x, y=0) = tf(x, y=z=0). Pions with vanishing transverse momentum (left column) are seen to come 



from a broad region symmetric around the center and are emitted rather late. Pions with K± = 0.5 GeV 
pointing in x-direction, on the other hand, are emitted on average somewhat earlier and only from a rather 
thin, crescent shaped sliver along the surface of the fireball at its point of largest transverse extension. 
The reason for this apparent "opacity" (surface dominated emission) of high-p^ particles is that they 
profit most from the radial collective flow which is largest near the fireball surface. Low-p± pions don't 
need the collective flow boost and are preferably emitted from smaller radii (where the flow velocity 
is smaller) when the freeze-out surface eventually reaches these points during the final stage of the 
decoupling process. 
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Fig. 21: Longitudinal cuts through the pion emission function S(x, K), integrated over transverse coordinates, for Y — pions 
with two different values of K± as indicated [58]. 

Fig. 21 displays the time structure of the emission process along the beam direction. Especially 
for low tansverse momenta one clearly sees a very long emission duration, as measured in the laboratory 
(center-of-momentum frame). The reason is that, according to the assumed longitudinal boost invariance, 
freeze-out happens at constant proper time r = t 2 — z 2 and extends over a significant range in longitudinal 
position z. This range is controlled by the competition between the longitudinal expansion velocity 
gradient (which makes emission of Y = pions from points at large z values unlikely) and the thermal 
velocity smearing. If freeze-out happens late (large r), the longitudinal velocity gradient ~ f/r is 
small and pions with zero longitudinal momentum are emitted with significant probability even from 
large values of \z\, i.e. very late in coordinate time t. Note that this is also visible in the lower left 
panel of Fig. 20 where significant particle emission still happens at times where the matter at z = 
has already fully decoupled. We will see shortly that this poses a problem when compared with the 
data. The long tails at large values of \z\ and t can only be avoided by reducing tj (thereby increasing 
the longitudinal velocity gradient and reducing the z -range which contributes Y = pions) and/or by 
additionally breaking longitudinal boost-invariance by reducing the particle density or postulating earlier 
freeze-out at larger space-time rapidities \r]\ [72]. 

From earlier hydrodynamic calculations [123] it was expected that a fireball evolving through the 
quark-hadron phase transition would emit pions over a long time period, resulting in a large contribution 
Pj_(t 2 ) to the outward HBT radius and a large ratio R a /R s . This should be a clear signal of the time- 
delay induced by the phase transition. It was therefore a big surprise when the first RHIC HBT data 
[124, 125] yielded R /R s ~ 1 in the entire accessible K± region (up to 0.7 GeV/c). In the meantime this 
finding has been shown to hold true out to K± > 1.2 GeV/c [126]. 

Figure 22 shows a comparison of the experimental data with results from hydrodynamic calcu- 
lations [64]. Clearly, a purely hydrodynamic description with default initial conditions (solid lines in 
Fig. 22) fails to describe the measured HBT radii. The longitudinal and outward radii Ri and R Q are 
too large, R s is too small, and the K± -dependence of both R a and R s are too weak in the model. Since 
for longitudinally boost-invariant sources Ri is entirely controlled by the longitudinal velocity gradient 




Fig. 22: HBT radii from hydrodynamic calculations [64] (solid lines) together with data from STAR [124] and PHENIX [125]. 
The dotted lines give hydrodynamic radii calculated directly after hadronization whereas the other two lines refer to different 
assumptions about the initial conditions (see text). 



at freeze-out which decreases as 1/t/, making Ri smaller within the hydrodynamic approach requires 
letting the fireball decouple earlier [64] or breaking the boost invariance [72]. As noted above, a smaller 
Ri would also reduce the emission duration, i.e. (t 2 ), and thus help to bring down R a , especially if 
freeze-out at nonzero rapidities \rj\ > happens earlier than at midrapidity. 

The fireball can be forced to freeze out earlier by changing the freeze-out condition (e.g. by im- 
posing freeze-out directly at hadronization, see dotted line in Fig. 22). This generates, however, serious 
conflicts with the single-particle spectra (see Sec. 4.1.9). Alternatively, one can allow transverse flow to 
build up sooner, either by seeding it with a non-zero value already at r cq (short dashed line labelled "hy- 
dro with FS" [64, 73]) or by letting the hydrodynamic stage begin even earlier (e.g. at Tf orm = 0.2 fm/c, 
long-dashed line). The last two options produce similar results, but do not fully resolve the problems 
with the magnitudes of Ri and R Q . 

The situation may improve by taking also the breaking of longitudinal boost-invariance into ac- 
count [72], but a fully consistent hydrodynamic description has not yet been found. In particular, the 
sideward radius R s is still too small and the Kx-depenences of both R s and R a are still to weak [64, 72]. 
It is hard to see how to increase R s without also increasing R Q and Ri which are already too large. Hy- 
brid calculations [127] in which the hydrodynamic Cooper-Frye freeze-out is replaced by transition to 
a hadronic cascade at T cr , followed by self-consistent kinetic freeze-out, tend to increase R s by making 
freeze-out more "fuzzy", but at the expense of also increasing R Q and Ri in a disproportionate manner, 
mainly due to an increase in the emission duration. This makes the problems with the R /R s ratio even 
worse. 

In the past, a strong Kj_ -dependence of R s has been associated with strong transverse flow [103, 
102]. It is therefore surprising that even the hydrodynamic model with its strong radial flow cannot 
reproduce the strong Kj_ -dependence of R s measured at RHIC. Also, according to Eqs. (55,56) the 
difference between R 2 , and R 2 S can be reduced, especially at large K± [121], if the positive contribution 
from the emission duration (i 2 ) is compensated by "source opacity", i.e. by a strongly surface-dominated 
emission process [115, 116, 128]. In this case the geometric contribution (x 2 ) to R 2 is much smaller 
than R 2 S = (x 2 ). Again, the source produced by the hydrodynamic model (see top right panel in Fig. 20) 
is about as "opaque" as one can imagine [51], and it will be difficult to further increase the difference 
(x 2 -x 2 ) [128]. 

This leaves almost only one way out of the "HBT puzzle", namely the space-time correlation term 



— 2/3_i_ (x t) in expression (56) for B%. It correlates the freeze-out position along the outward direction 
with the freeze-out time. The hydrodynamic model has the generic feature that, in the region where most 
particles are emitted (see Fig. 20), these two quantities are negatively correlated, because the freeze-out 
surface moves from the outside towards the center rather than the other way around. Hence the term 
— 2/3j_ (x t) is positive and tends to make i? 2 , larger than R 2 S . The small measured ratio R /R s < 1 may 
instead call for strong positive x a —t correlations, implying that particles emitted from larger x Q values 
decouple later. Hydrodynamics can not produce such a positive x Q —t correlation (at least not at RHIC 
energies). On the other hand, there are indications that microscopic models, such as the AMPT [122] 
and MPC [129] models, may produce them, for reasons which are not yet completely understood. 

One should also not forget that Fig. 22 really compares two different things: The data are extracted 
from the width of the 2-particle correlator in momentum space while in theory one calculates the same 
quantities from the source width parameters in coordinate space. The two produce identical results only 
for Gaussian sources. The hydrodynamic source function shown in Figs. 20 and 21 are not very good 
Gaussians and show a lot of additional structure. We have checked, however, by explicit computation of 
the momentum-space correlation function using Eq. (49) that the non-Gaussian effects are small. 2 The 
largest non-Gaussian effects are seen in the longitudinal radius Ri [112], but although the corresponding 
corrections go in the right direction by making the Ri extracted from the momentum-space correlator 
smaller, the effect is only a fraction of 1 fm and not large enough to bridge the discrepancy with the data. 

It was recently suggested [94, 130] that neglecting dissipative effects might be at the origin of 
the discrepancy between the purely hydrodynamic calculations and the data. A calculation of first-order 
dissipative corrections to the spectra and HBT radii at freeze-out [94], with a "reasonable" value for the 
viscosity, yielded a significant decrease of Ri along with a corresponding strong reduction of the emission 
duration contribution to R Q , both as desired by the data. There was no effect on the x D —t correlations, 
however, and only a weak effect on R s which went in the wrong direction, making it even flatter as a 
function of K±. The rather steep K± -dependence of the data for both R s and R Q and the larger than 
predicted size of R s at low K± are therefore not explained by this mechanism [94, 129]. Furthermore, the 
elliptic flow V2 has been shown to be very sensitive to viscosity [94, 95], and the viscosity values needed 
to produce the desired reduction in Ri turned out [94] to reduce V2 almost by a factor 2, incompatible 
with the data. The "RHIC HBT puzzle" thus still awaits its resolution. 

5.5 Non-central (b^O) collisions — azimuthally sensitive HBT interferometry (asHBT) 

In non-central collisions the azimuthal symmetry of the source is broken by the impact parameter vector 
b which, together with the beam direction, defines the reaction plane. Pairs emitted at different azimuthal 
angles $ = Z(6, K±) relative to the reaction plane therefore originate in general from effective emission 
regions (homogeneity regions) of different sizes and shapes. 

As a result, the components of the spatial correlation tensor S^ u = {x^Xy) now depend not only 
on the magnitude of the transverse pair momentum K±, but also on its angle $ relative to the reactuion 
plane. Furthermore, the source is no longer reflection symmetric with respect to the plane defined by 
the beam axis and K± (i.e. symmetric under x s — ► — x s ), so R os and i?; s no longer vanish. This means 
that for non-central heavy-ion collisions we have to measure 6 (not 4) HBT radii, and they are functions 
not only of Y and K±, but also of the emission angle relative to the reaction plane. A complete HBT 
analysis of non-central heavy-ion collisions thus requires the event-wise determination of the reaction 
plane (through a Fourier analysis of the single particle spectrum as, for example, used for extracting the 
elliptic flow) and the investigation of the dependence of the HBT radii on the angle <I> relative to that 
plane. 

2 They would have been larger if we had included resonance decay pions in the emission function, as found by Lin et al. 
[122], which was our main reason for not doing so. 



5.5.1 asHBT formalism 

Experimentally it is still useful to characterize the source in terms of HBT radii R Q and R s parallel and 
perpendicular to the transverse emission direction defined by the pair momentum K±. The emitting 
source S(x, K), however, is most conveniently written down in transverse coordinates x, y aligned with 
the reaction plane, since it typically (at least for collisions between spherical nuclei) is symmetric under 
reflection with respect to that plane, y^—y combined with — ► — The analogue relations to Eqs. (55- 
57), which express the HBT radii through the components = {x^x u ), therefore involve an explicit 
rotation by the angle between the (x,y) axes characterizing the source S^ v and the (x Q ,x s ) axes 
defined by the emission direction [131, 132, 133]: 

R 2 S = \{S XX + Syy) ~ \{S XX Syy) CO S (2$) - S X y S ili{2^) 
R 2 = \{S XX + Syy) + \{S XX -Syy) COS(2$) + S X y fflll(2$) 

-2(3 ± (S tx cos $+S ty sin $) + p 2 ± S tt , 
Rls = S xy cos(2$) - \ (S xx -S yy ) sin(2$) + ± (S tx sin <f>-S ty cos $), 
R? = S zz - 2(3 L S tz + PiSu, 

Rli = {S xz -p L S tx ) cos® + (S yz -(3 L S t y) sin® -f3 ± S tz + p L P±Stt, 

R 2 sl = (S yz -(3 L S ty ) cos® -(S xz -(3 L S tx ) sin®. (75) 

These equations exhibit the explicit $ -dependence resulting from the rotation between the (x, y) and 
(x ,x s ) axes, but not the implicit one reflecting the $ dependence of the spatial correlation tensor 
S tlu (Y, K±, The total emission angle dependence of the HBT radii results from the combination 
of both. 

The implicit ^-dependence of the spatial correlation tensor is restricted by symmetries of the 
source [133]. It is a relativistic effect associated with an azimuthal spatial source deformation superim- 
posed by strong transverse collective flow [131, 134] which vanishes with the 4 th power of the transverse 
flow velocity v±/c for weak or no collective expansion [132, 134]. 

A full analysis of the symmetry constraints on S I1U (Y, K±, <£) for symmetric collisions between 
spherical nuclei and for pairs detected in a symmetric rapidity window around Y = was performed in 
Ref. [133]. One finds the following most general form for the azimuthal oscillations of the HBT radii: 
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We see that only even or odd sine or cosine terms occur, but no mixtures of such terms. Statistical er- 
rors in the resolution of the reaction plane angle as well as finite angular bin sizes in <E> tend to reduce 
the actually measured oscillation amplitudes; fortunately, these dilution effects can be fully corrected 
by a model-independent correction algorithm [133, 135]. A Gaussian fit to the thus corrected correla- 
tion function, binned in Y, K± and emission angle <!>, then yields the "true" HBT radius parameters 



Rfj(Y, K±, <J>) from which the n th order azimuthal oscillation amplitudes are extracted via 



i "bin 

4- n (Y, K ± ) = J2 4( y > K ^ ^)osc(n^). (77) 

^bin j = i 

Here n-bin indicates the number of (equally spaced) $ bins in the data and osc(n$j) stands for sin(n<I>j) 
or cos(n$j) as appropriate, see Eqs. (76). (Note that Nyquist's theorem limits the number of harmonics 
that can be extracted to n < n^ n .) 

We would like to relate the azimuthal oscillation amplitudes of the 6 HBT radius parameters to 
the geometric and dynamical anisotropics of the source, as reflected in the azimuthal oscillations of 
the 10 independent components of the spatial correlation tensor. Their allowed oscillation patterns at 
midrapidity Y = are given by [133] 

A($) = \{x 2 +y 2 ) = A)+2En>2,cven^nCO S (n$), 

= \(x 2 -y 2 ) = 5o+2E„> 2 ,eve„5„co S (n$), 
C($) = (xy) = 2£ re >2,evenCr i sin(n$), 

= = L> +2E„>2,even^nCOs(n$), 

E($) = {tx) = 2E„>i,odd^nCos(n*), 
F($) = (ty) = 2E„>i,oddi ? nSin(n$), 
G(*) = {tz) = 2E„>i,oddG , nCos(n$), 

H($) = {XZ) = iJ +2En>2,even^nCOs(n$), 
/($) = (P) = 2E„>2,even-^nCOs(n$), 

J($)EE(5 2 ) = J +2E n >2,even^nCOs(n$). (78) 

The missing terms in the sums over n have amplitudes which are odd functions of Y and vanish at midra- 
pidity. They do, however, contribute to the HBT radii if the data are averaged over a finite, symmetric 
rapidity window around Y = [133]. Their contributions can be eliminated by varying the width AY of 
this rapidity window and extrapolating quadratically to Ay — ► 0. Note that Cq = Eq = Fq = Gq = Iq = 
by symmetry, i.e. the corresponding components of oscillate around zero. 

The oscillation amplitudes of the HBT radii relate to the oscillation amplitudes of the source 
parameters as follows [133, 136]: For the odd harmonics n= 1, 3, 5, ... we have 
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Rsl,n = ^(—Hn-l+H n +l+In-l+In+l — PL(—E n -i+E n+ i+F n -i+F n+ i)}, 

whereas the even harmonics n = 0, 2, 4, . . . satisfy 

Rs,n = (An) + \(-F> n ~2-B n+ 2+C n -2 — C n+ 2), 
^o,n = (Ai) + \(Bn-2+B n+ 2-C n -2+C n+ 2) 

—(3±(E n -i+E n+ i—F n -i+F n+ i) + (3j_(D n ), 

F^>s,n = ^(—Bn-2+B n+ 2+C n -2+Cn+2+P±(E n -l—E n+ i — F n -i—F n+ i)), 

R\ n = (Jn)-2(P L Gn) + (P 2 L D n ). (80) 

In these relations it is understood that all negative harmonic coefficients n < as well as Cq,Eq,Fq, Go 
and Iq are zero. The angular brackets (...) indicate an average over a finite, symmetric rapidity window 
around Y = 0. The terms involving the longitudinal pair velocity [3l vanish quadratically as the width 
Ay of that window shrinks to zero. 



Even after extrapolating to Y = in this way, we have still many more source parameters than 
measurable HBT amplitudes. One counts easily that up to n = 2 there are 9 measurable Fourier coef- 
ficients which (at Y = 0) depend on 19 source amplitudes. From there on, increasing n by 2 yields 6 
additional measured amplitudes which depend on 10 additional source amplitudes. This lack of analysis 
power is an intrinsic weakness of the HBT microscope and due to the fundamental restrictions arising 
from the mass-shell constraint q° = (3 ■ q. The reconstruction of the source thus must necessarily rely on 
additional assumptions. 

One such assumption which may not be too unreasonable is that the emission duration D = (t 2 ) is 
approximately independent of emission angle and that the source is sufficiently smooth that higher order 
harmonics n > 3 of S^ v can be neglected. Such source properties would result in the "Wiedemann sum 
rule" [131] 

^,2-^,2 + 2^,2 = (81) 

which can be experimentally tested. If verified for all K± it would provide strong support for the under- 
lying assumptions on the source. In this case we can measure 3 azimuthally averaged HBT radii and 5 
independent oscillation amplitudes with n < 2, depending on 14 source parameters of which 5 can be 
eliminated by going to K± = (3± = (see [133] for explicit expressions). This makes the geometry of the 
effective source for particles with K± = "almost solvable" 3 , as confirmed by hydrodynamical calcula- 
tions [5 1 ] which show that at K± = the effective emission region closely tracks the overall geometry 
of the source even if it is strongly and anisotropically expanding. 

The source geometry can be completely reconstructed from HBT data if transverse flow is so weak 
that all implicit $ -dependence (i.e. all higher harmonics n > 1) of can be neglected. In this case 
one obtains at Y = the "geometric relations" [132] 

R 2 fi = A = ±(x 2 +f) , 

R l,o ~ r2 s,o = P± D o = /3±(P)o, 

R 2 ,o = 'h = (z 2 )o, 

R ll,l = ~ r2 sI,1 = k H ° = l(^)0; 

Rl,2 = ~R 2 s,2 = ~R 2 os,2 = \B = l(x 2 -y 2 ) . (82) 

Ao describes the average transverse size and Bq (which generates a second-order harmic in the transverse 
HBT radii) the transverse deformation of the source. Hq generates a first-order harmonic in the ol and si 
cross terms and describes a longitudinal tilt of the source away from the beam direction [132]. Such a tilt 
was found in Au+Au collisions at the AGS [138]. Its sign yielded important information on the kinetic 
pion production mechanism [132, 138]. 

5.5.2 Azimuthal oscillations of HBT radii in hydrodynamics and RHIC data 

We can use the hydrodynamic model introduced in Section 4.1.9 to describe the momentum spectra 
and explore what it would predict for the azimuthal oscillations of the HBT radii in non-central Au+Au 
collisions at RHIC [51]. The most interesting question, perhaps, is whether we can reconstruct the 
orientation and magnitude of the spatial deformation of the source at freeze-out from the HBT analysis, 
thereby complementing our knowledge of the orientation and magnitude of the momentum anisotropy 
from the elliptic flow analysis. 

3 Note that the limit K± — » eliminates all influence from the temporal structure of the source; the emission duration and 
correlations between position and time at freeze-out must be extracted from correlation data at non-zero K± where, however, 
the explicit (3± -dependence associated with factors t in the variances ix^x-J) interferes with the implicit K± -dependence in 
the spatial variances which result from collective expansion flow [103]. Models for disentangling these different contributions 
to the K± -dependence have been extensively studied for azimuthally symmetric sources [103], but must be generalized for 
azimuthally deformed sources [137]. 



To see how different orientations of the spatial deformation (in-plane elongated vs. out-of-plane 
elongated) manifest themselves in the azimuthal oscillations of the HBT radii, one can explicitly con- 
struct sources of either orientation and perform an HBT analysis on them. Since the elliptic flow makes 
the fireball expand faster into the reaction plane than perpendiculer to it, we can generate an in-plane 
elongated source (IPES) at freeze-out from the initial out-of-plane oriented overlap region simply by 
making sure that the elliptic flow has enough time to act before the system decouples. By increasing the 
initial energy density or temperature, we can make the total time until decoupling arbitrarily large. Let 
us therefore study hydrodynamic sources evolving from different sets of initial conditions [51]. 

The first set (labelled "RHIC1") corresponds to optimized values for central Au+Au collisions at 
y / SNN = 130 GeV which were used earlier in Section 4.1.9 (see Table 1). The second set of initial condi- 
tions (labelled "IPES") uses a much higher initial temperature T eq = 2 GeV (at b = 0), and the hydrody- 
namic evolution is started at r eq = 0.1 fm/c and stopped when a freeze-out temperature Td ec = 100 MeV 
has been reached. Such a high initial temperature can probably not even be achieved at the LHC, and 
the strong transverse flow generated in this case probably causes the system to decouple already much 
closer to the hadronization temperature T cr = 164 MeV. But if we force the system to start at such a high 
temperature and freeze out so low, we give it enough time to convert the initial out-of-plane deformation 
into a final in-plane deformation (which is the phenomenon we want to study). The following analysis 
refers to semi -peripheral Au+Au collisions with impact parameter 6 = 7 fm. 
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Fig. 23: (a) Cuts through the freeze-out surface tf (r) at z — along and perpendicular to the reaction plane, (b) Contour plots 
in the transverse plane of the time-, z-, and momentum-integrated emission function d 2 N/d 2 r for RHIC1 and IPES initial 
conditions (see text), (c) Cuts through diagram (b) along and perpendicular to the reaction plane, for RHIC1 and IPES initial 
conditions. In all cases 6 = 7 fm. 



Figure 23a shows the freeze-out surface tdec(^) m cut s along (y = 0) and perpendicular to the 




Fig. 24: Contours of constant emission density in the transverse plane for RHIC1 initial conditions (see text). The thick line 
indicates the largest transverse extension of the freeze-out hypersurface (see Figs. 23a,b). The four panels show emission 
regions for midrapidity pions (Y — 0) with K± = and, for three emission angles indicated by arrows, with K± — 0.5 GeV. 

reaction plane (x = 0). Initially the source is extended out-of-plane (larger in y than in x direction), but it 
then expands more rapidly into the x-direction, becoming in-plane elongated at later times. For RHIC1 
initial conditions this only happens after most of the matter has already decoupled; as a consequence 
the time-integrated source, shown in Fig. 23b, is still longer in y than in x direction in the RHIC1 case. 
For IPES initial conditions the deformation changes sign before most particles decouple, and the time- 
integrated source appears in-plane-extended (see again Fig. 23b). Also, it is much larger due to the much 
higher initial energy density and longer lifetime. 

Figure 23c shows cuts along and perpendicular to the reaction plane through the density contour 
plots d 2 N/d 2 r = J (d 3 K/Ex)dz dt S(x, K) of Fig. 23b. Pion emission is seen to be strongly surface 
peaked, in particular at RHIC1 where the freeze-out radius is almost constant for a long time. This "opac- 
ity" is weaker both at lower collision energies (where the freeze-out surface shrinks to zero continuously 
[50]) and at higher energies, due to larger temporal variations of the freeze-out radius. 

Figures 24 and 25 show the spatial distributions of pions emitted with fixed momentum. Shown 
are density contours of / dz dt S(x, K) in the transverse plane (x, y) for pions with rapidity Y = and 
fixed K±, for three emission angles $ = 0, 45° and 90° relative to the reaction plane. (In Fig. 25 we 
replaced <I> = 45° for clarity by the equivalent angle = 225°.) Particles with vanishing transverse mo- 
mentum K± are seen to be emitted from almost the entire interior of the "bathtub" shown in Figs. 23b,c; 
for RHIC1 (IPES) initial conditions this region is elongated out-of-plane (in-plane). For slow pions the 
source thus looks transparent. Pions with sufficiently large transverse momenta are emitted from rela- 
tively thin regions close to the rim of the "bathtub" where the flow velocity is largest and points into the 
direction of the emitted pions. For fast pions the source thus looks opaque, more so at RHIC than at 




Fig. 25: Same as Fig. 24, but for IPES initial conditions. 

higher energies. Their emission regions rotate with the emission angle, constrained by the shape of the 
"bathtub". 

The HBT radii, calculated from the widths of the i^-dependent emission regions according to 
Eqs. (75), are shown in Fig. 26, as functions of the azimuthal emission angle <P. (Since Y = 0, all 
terms ~ vanish.) Note that R 2 Q and R 2 S receive purely geometric and mixed space-time correlation 
contributions; the latter are proportional to the pair velocity (3±. Since we want to use pion interferometry 
to obtain information on the geometric deformation of the source at freeze-out, these should be analyzed 
separately. In the right panel of Fig. 26 we therefore show the geometric contributions separately as thin 
circled lines. For a detailed discussion of all the curves in Fig. 26 I refer the reader to Ref. [51]; here I 
will only discuss the most important point. 
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Fig. 26: Left: Azimuthal oscillations of the HBT radii at Y = for b = 7 fm Au+Au collisions at s/sim = 130 GeV (RHIC 1), 
for four values of the transverse momentum K± as indicated. Right: The same for IPES initial conditions. For R% and R^, s the 
geometric contributions are shown separately as thin circled lines. 

Let us first look at the geometric contributions. From Figs. 24 and 25 we already know that pions 
with K± = are emitted from almost the entire fireball and thus probe the different sign of the spatial 
deformation of the total (momentum-integrated) RHIC1 and IPES sources shown in Fig. 23. This is 
reflected by the opposite sign of the oscillation amplitudes of R 2 S and of the geometric contribution to 
R 2 at K± = in the left and right panels of Fig. 26. At higher K± -values the oscillations for R 2 S in the 
right panel Fig. 26 change sign, but those of the geometric contribution to R 2 do not. This reflects an 



intricate interplay between geometric and flow effects [51]. However, the opposite signs of the oscillation 
amplitudes near K± = of R 2 , R 2 , and R s 2 for the RHIC1 and IPES sources are clear and unique 
indicators of the opposite directions of the spatial deformation for these two sources and can thus be 
used to diagnose the latter. In particular, the ratio of the second harmonic oscillation amplitude R 2 S 2 of 
the sideward radius to its <£-averaged value R 2 is seen to be almost independent of K± [51, 137]. This 
ratio can therefore be used to extract the spatial deformation of the entire source even at non-zero (albeit 
not too large) K±, even though the homogeneity regions at non-zero K± do not probe the entire source 
[137, 139]. 
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Fig. 27: Left: Azimuthal oscillations of the HBT radii near Y = from Au+Au collisions at = 200 GeV, for 20%-30% 

centrality events and four bins of the transverse momentum K± (GeV/c) [139]. The lines show fits with Eqs. (76). Right: 
Zeroth order (left column) and normalized second order (right column) harmonic coefficients of the HBT radii, for three K± 
bins (stars, circles, and triangles), as a function of collision centrality, from the same STAR experiment [139]. Larger participant 
numbers correspond to more central collisions. 

Recently, the STAR Collaboration published the first azimuthally sensitive HBT analysis of Au+Au 
collisions at RHIC [139]. As you can see in Fig. 27, the quality and level of detail of the data is impres- 
sive. The analysis is fully 3-dimensional in the relative momentum q and sufficiently finely binned in 
K±, <£, and collision centrality to obtain a reasonably complete HBT snapshot of the final state of the 
collision fireball created at RHIC. (Of course, theorists will always keep asking for more detail, but even 
that is coming.) It is clear that with these asHBT data a new dimension has opened, with unprecedented 
constraining power on theories ! 

Comparing the left panels in Figs. 26 and 27 one sees very strong similarity. The absolute scales on 
the vertical axes don't match, but this just reflects the "RHIC HBT puzzle" from the previous subsection, 
namely that the hydrodynamic model does not reproduce the angle-averaged HBT radii and overpredicts 
Ri and R a while underpredicting R s . The signs of the oscillation amplitudes match, however, clearly 
indicating that the measured source is still somewhat out-of-plane elongated at freeze-out, as predicted 
by the hydrodynamic model (see Figs. 23 and 24). The centrality systematics shown in the right panel of 
Fig. 27 shows that the HBT radii grow as the collisions become more central (left column - no surprise 
there), but that the growth rate decreases, together with the radii themselves, as the pair momentum K± 
increases. This is qualitatively consistent with strong longitudinal and radial flow, as discussed in the 



previous section, even though the measured K± -dependence of R 2 and R 2 is stronger than predicted. 
The experimental data show that the strong K± -dependence of R 2 in particular is not only a feature of 
central collisions, but seems to persist towards peripheral collisions. 

The rightmost column of Fig. 27 shows that the oscillation amplitudes are are largest in peripheral 
collisions and decrease to zero in central collisions. Again this is no surprise since the oscillations reflect 
the geometric deformation of the source and peripheral collisions produce, at least initially, the most 
deformed fireballs. Note the very weak K± -dependence of the relative oscillation amplitudes which is, at 
least qualitatively, consistent with hydrodynamic predictions. A more quantitative analysis of the source 
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Fig. 28: Final source eccentricity at freeze-out, £fi na i, extracted from the normalized second order harmonic coefficient of 
the sideward HBT radius, Rs,2/Rl,o [137] (see text), as a function of the initial eccentricity, einitiai, from a Glauber model. 
Smaller eccentricities correspond to more central collisions. The dotted line indicates Eimtiai = £fi na i- The data are from Au+Au 
collisions at ^/s^= 200 GeV meassured by STAR [139]. 

deformation as a function of centrality is presented in Fig. 28 from the same publication [139], where 
the final spatial eccentricity £fi na i at decoupling, extracted from the normalized oscillation amplitude 
R 2 2 /R 2 with the help of a blast- wave model parametrization of the source [137], is contrasted with the 
initial spatial deformation £i n iti a i of the original nuclear overlap zone, calculated from the Glauber model. 
At all impact parameters, the final deformation is smaller than the initial one, showing the consequences 
of elliptic flow. Relatively speaking, the reduction of the eccentricity is stronger in almost central than 
in peripheral collisions, but at no centrality does the freeze-out configuration exhibit a negative value 
for £fi na i, which would indicate an in-plane-elongated source. For mid-peripheral (b ~ 7fm) collisions 
(second point from the right), the extracted value £fi na i = 0.11 ± 0.035 agrees (within the relatively large 
systematic error of about 30%) with the eccentricity of the hydrodynamic source shown in Figs. 23,24 
which has e x ~ 0.14. This provides further support for the hydrodynamic model, since it indicates an 
intrinsic dynamical consistency between the evolution of the eccentricity from its initial to its final value 
and the hydrodynamic expansion time scales and flow profiles. A stronger and more confident statement 
requires, of course, the resolution of the "HBT puzzle", i.e. a quantitative theoretical understanding of 
the angle-averaged HBT radii which is so far missing. 

5.6 What has HBT interferometry of heavy-ion collisions taught us? 

Let me quickly summarize this long chapter by stating the following as our Third Lesson: Two-particle 
momentum correlations, in particular Bose-Einstein correlations between pairs of identical bosons, are 
a powerful tool to constrain the space-time structure and evolution of the heavy-ion fireball, thereby 
yielding complementary information to its momentum-space structure extracted from the single-particle 
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spectra. Emission-angle dependent HBT interferometry for non-central collisions (or central collisions 
between deformed nuclei) allows to extract the shape and orientation of the deformed source at freeze- 
out. By comparing it with the initial spatial deformation of the nuclear overlap region when the nuclei 
first hit each other, we get an idea how the system evolved and can constrain the time between impact 
and decoupling. Azimuthally sensitive HBT interferometry (asHBT) is a newly developing field, but 
first results indicate consistency of the extracted final deformation with hydrodynamic evolution mod- 
els. On the other hand, the absolute magnitudes of the emission-angle averaged HBT radii and their 
dependence on the transverse momentum K± of the emitted pairs is still poorly understood: Almost all 
dynamical models overpredict Ri and R Q , underpredict R s and are unable to reproduce the measured 
strong K± -dependence of the transverse radii R s and R Q . This "RHIC HBT Puzzle" indicates that our 
understanding of the kinetics of the freeze-out process is still insufficient and constitutes, in my opinion, 
the most important open problem in our present global picture of heavy-ion collision dynamics. 

6 STATISTICAL HADRONIZATION AND PRIMORDIAL HADROSYNTHESIS 

In this section we'll discuss what kind of information can be extracted from the measured hadron abun- 
dances (i.e. from the normalization of the hadron spectra, in contrast to their shape which we discussed in 
the previous section). Up to this point we did not discuss the normalization of any of the spectra, nor did 
we care about it. All we needed for local thermal equilibrium and the build-up of thermal pressure and 
collective flow was rapid momentum equilibration, and this can in principle be achieved entirely by elas- 
tic collisions (including resonant elastic collisions in the hadronic phase such as ir+N — > A — > ir+N) 
which do not change any of the particle abundances. Whether or not the microscopic interactions alter 
the chemical composition or leave it unchanged is, for this part of the discussion, entirely irrelevant. 

This analysis of the spectral shapes has already taught us that hadrons are not formed right away 
but appear with some delay via the hadronization of a dense, thermalized, collectively expanding quark- 
gluon system. Let us now consider this hadronization process and ask how it affects the hadron abun- 
dances, i.e. the chemical composition of the hadronic system emerging from this process. 

6.1 Statistical hadronization and the Maximum Entropy Principle 

The microscopic description of the process in which thousand of quarks and gluons combine to form 
thousands of final state hadrons is clearly an impractical problem. But it is also not needed. A micro- 
scopic approach is asked for if we want to understand the (relatively rare) processes by which, say, a 
high-p^ hadron of a given flavor is created. Such particles probe the microscopic dynamics of high-p^ 
partem production, which is controlled by QCD, and although they may rescatter in the dense fireball 
medium they never completely lose their memory of the primary production process. We can exploit this 
and use them as "deep probes" for the fireball matter, but extracting the information on their production 
and subsequent medium modification then clearly requires a detailed microscopic treatment, including 
the hadronization of such hard partons. This is completely different for soft partons and hadrons: in their 
case the rescattering by the medium wipes out all memories about their initial production process, and 
there is no way their final momenta could ever be traced back to their initial momenta. This is in partic- 
ular true for the hadronization process which is probably the most non-perturbative process imaginable: 
there are literally innumerable different microscopic channels which might lead to the production of, say, 
a negative pion with a soft momentum of order 300 MeV/c in the laboratory. 

The proper approach to such a non-perturbative, multi-channel multi-particle problem is (as al- 
ready recognized by Fermi, Landau and Hagedorn more than half a century ago) a statistical one. From 
a microscopic point of view, there is really very little information we can extract from the yields and 
momentum spectra of soft single hadrons. Due to the strong interactions among the quarks and gluons as 
they are forming hadrons, anything that is not explicitly forbidden by the conservation laws for energy, 
momentum, baryon number, electric charge, and net strangeness will in fact happen. (Net strangeness can 



only be changed by weak interactions, on a time scale which is much too long compared to a heavy-ion 
collision.) 

Note that such a statistical approach has, of course, its limitations: when we start to investigate 
correlations between pairs, triplets, quadruplets etc. of particles, at some level they will begin to realize 
that they all belong to a single heavy-ion collision event and are not produced entirely inependently, 
but correlated in a non-statistical way. (For example, the momenta of all 2479 particles emerging from 
the collision have to add up to the total momentum of the initially colliding nuclei, and this is likely 
to generate non-statistical momentum correlations among considerably smaller subclusters of particles.) 
The statistical approach is certainly best justified for single-particle observables, such as hadron spectra 
and abundances. 

These words are translated into equations via the Maximum Entropy Principle [140]: Let us con- 
sider the entropy S of a small fluid cell of size AV in its own rest system, 



d 3 x d?p 
(2vr) 3 



Si In Si + <?i(l - In (1 - ViSi) , (83) 



where the sum over % includes all particle species in the cell and i?j = + 1(— 1) for fermions (bosons). We 
assume that soft hadronization happens through many different channels which are constrained only by 
the local (in AV) conservation laws for energy, baryon number and strangeness (J doj = J AV J ^j0?r- )'- 

(E)= [cLu^EiSi, (B)= [<koJ2BiSi, (S) = JdLoY,Sifr = 0. (84) 

i i ' i 

Here Bi and S% are the baryon number and strangeness carried by each particle of species i, and we used 
that the net strangeness of the two colliding nuclei is zero and that strong interactions keep it that way 
locally. Maximizing S with these constraints gives for the most likely distribution for particle species i 
[141] 

f l (E l ,AV) = emv)[E ^ MAV)]+ ^ with Vi(AV) = B i » B (AV) + S i » s (AV), (85) 

where (3, hb and fis are Lagrange multipliers related to the values of the constraints (84) in Ay. If 
the hadronization process does not change the number of ss pairs (for example, because strange valence 
quarks are too heavy), we should implement an additional constraint on the total number of strange 
quarks and antiquarks which fixes that number to the value before hadronization begins. In this case 
the expression for /ij in (85) generalizes to //j = BifiB+SiHs+\si\p, 8 where is the total number of 
strange quarks and antiquarks in hadron i and jl s is a further Lagrange multiplier related to this additional 
constraint. The associated fugacity 7 S = e^ s is known as the "strangeness saturation factor" [142]. An 
under- resp. oversaturation of the strange particle phase space relative to its chemical equilibrium value 
corresponds to 7 S < 1 resp. 7 S > 1. 

Equation (85) is a local thermal and chemical equilibrium distribution function. If the system could 
be kept at constant volume, any type of strong interaction among the hadrons would leave this distribution 
unchanged since such microscopic processes again conserve energy, baryon number and strangeness. 
However, the apparent "equilibrium" expressed by Eq. (85) is not achieved kinetically (i.e. as a result of 
hadronic rescattering) but statistically (by interference of many different hadron production channels). It 
does not require a thermalized prehadronic state (although very strong deviations from equilibrium, e.g. 
at high p±, in the prehadronic state may survive the hadronization process and invalidate the Maximum 
Entropy Principle). 

It is very important to realize that the same local equilibrium distribution (85) can be the result 
of two conceptually entirely different types of processes: we can either take system of hadrons with an 
arbitrary (except for the constraint on total energy, baryon number and strangeness) initial phase-space 



distribution and let it evolve for a sufficiently long time to obtain Eq. (85) as a result of the action of 
elastic and inelastic processes among the hadrons. This is kinetic equilibration. Or we produce the 
system of hadrons of given energy, baryon number an strangeness from some non-hadronic state by a 
statistical process which fills hadronic phase-space in the statistically most probable configuration. This 
is statistical equilibrium. Both processes share the property that they lead to a state of Maximum Entropy. 
However, statistical hadronization can produce a Maximum Entropy distribution through non-hadronic 
processes which occur much faster than any inelastic scattering among hadrons at the given energy and 
baryon number density. 

If you try to create, via statistical hadronization, a "pre-established equilibrium" distribution (85) 
of hadrons at energy density e > e cr (T > T cr ), lattice QCD tells us that such a "superheated" hadronic 
equilibrium configuration is unstable against deconfinement. If you leave this system to itself, the 
hadrons will dissolve again into quarks and gluons. Statistical hadronization thus can only proceed once 
the energy density has dropped to ehad = £cr — 1 GeV/fm 3 . Therefore the hadrons, when first formed 
from the hadronizing quark-gluon state, will have apparent chemical equilibrium abundances correspond- 
ing to a temperature Th a d = T CT . A Before dicussing what happens to this "chemical temperature" after 
hadronization has been completed, let us first have look at the data. 



6.2 Hadrosynthesis at SPS and RHIC and pre-established hadronic chemical equilibrium 

In Figure 29 I show the particle multiplicity ratios for stable (with respect to strong interactions) hadron 
emitted in Pb+Pb collisions at the SPS (left) and in Au+Au collisions at RHIC (right). In both cases the 
data have been fit to an "apparent chemical equilibrium distribution" (85); the resulting fit parameters 
are given in the Figure. The left panel presents the quality of the fit in a form introduced by Becattini 
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Fig. 29: Chemical equilibrium fits to the hadron abundances from Pb+Pb collisions at -^/snn = 17GeV [54, 145] (left) and 
from Au+Au collisions at ^/snn = 17 GeV [68, 146]. The baryon chemical potential \ib used on the left and the quark 
chemical potential \i q used on the right are related by fiB = 3/i q . The strangeness chemical potential fi s is not an independent 
fit parameter, but follows from overall strangeness neutrality. The strangeness saturation factor 7„ [142] is approximately 1 in 
both cases. 

[144] where the measured multiplicities on the vertical axis are plotted against the multiplicities from 
the thermal model after optimization of the thermal parameters. If Eq. (85) gave a perfect description 
of the data, all points would lie exactly on the diagonal. For clarity the differences between data and 

4 I am not stating the same for the slope of their momentum distribution because the hadronizing quark-gluon matter already 
undergoes transverse flow which blueshifts the momentum spectrum. The abundances and the "chemical temperature" extracted 
from them are not affected by flow [143]. 



optimized thermal model are plotted (in terms of standard deviations) on a linear scale at the bottom. In 
the right panel the data are compared with bands of predicted ranges whose width is controlled by the 
statistical and systematic error of the thermal fit parameters [147]. Note that in the fit of the RHIC data 
only a fraction of the data was used for the fit while the rest (including ratios involving the rare doubly 
and triply strange S and hyperons) were successfully predicted by this fit. 

These fits tell us several things. First, not surprisingly, the baryon chemical potential at RHIC is a 
lot smaller than at the SPS, due to increased transparency of the target nucleus for the incoming projectile 
baryons. Second, the value for the strangeness saturation factor 7 S is almost equal to 1 in both cases; I 
will return to this observation in the following subsection. Third, and most importantly, the "chemical 
decoupling temperature" T c v lcm is in both cases much higher than the kinetic freeze-out temperature 
Td C c extracted from the shape of the momentum spectra. At the SPS T c h em ~ 160 MeV and at RHIC 
T c hem ~ 175 MeV. Interpreting r c h em i n the Maximum Entropy spirit as a Lagrange multiplier which tells 
us at which energy density statistical hadronization happened we see that the data corrspond to energy 
densities which are very, very close to the critical energy density for deconfinement predicted by lattice 
QCD. According to lattice QCD, hadrons in thermodynamic equilibrium cannot exist at temperatures 
higher than about 175 MeV if jjlb = 0, and this limit goes down if hb ^ 0. So these chemical decoupling 
temperatures are about as high as they could possibly be for a system of hadrons! 

Furthermore, we know from the spectra and elliptic flow that the fireball is collectively expanding 
and therefore undergoes rapid cooling. So we face the following conundrum: hadron formation is impos- 
sible at energy densities e corresponding to T > T c h cm > but once the hadrons are formed at e cr the system 
does not remain at this energy density but keeps diluting very rapidly. So if the hadrons first appeared 
with arbitrary abundances they would have had no time at all to adjust their abundances to chemical 
equilibrium values at T c h em ~ T CI by inelastic rescattering ! The possibility that the measured chemical 
freeze-out temperature T c h em — 160 — 170 MeV reflects a kinetically established chemical equilibrium is 
therefore logically excluded. The apparent chemical equilibrium at T c h em ~ T CT seen in the hadron abun- 
dances must thus be a "pre-established equilibrium" in the sense of statistical hadronization according to 
the Maximum Entropy Principle [15, 148]. 
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Fig. 30: (a) Pion multiplicity dN w /dy per unit overlap area S (stars), chemical decoupling temperature T ch (circles), thermal 
(kinetic) decoupling temperature Tki n (open triangles), and (b) the average transverse flow velocity (/3r) at kinetic decoupling 
as a function of collision centrality (measured by the total charged multiplicity pseudorapidity density dN cb dr/ at midrapidity) 
from 200 A GeV Au+Au collisions at RHIC [150], 

In addition to being logically excluded, there is also strong empirical evidence against a kinetic 
interpretation of the observed chemical decoupling temperature: First, it takes the same universal value 
of about 170 MeV in high energy e+e~ [144], pp and pp [15], and nucleus-nucleus collisions [149]. 



Second, in Au+Au collisions at RHIC, T c h cm is found to be independent of collision centrality (open 
circles in the left panel of Fig. 30), in contrast to the observed significant centrality dependence of the 
thermal decoupling temperature and average radial flow velocity extracted from the momentum spectra 
(triangles in Fig. 30 - see Fig. 9 for comparison and the discussion in Section 4.1.8). As discussed earlier, 
any kinetic decoupling process is a competition between local thermalization and global expansion rates. 
The latter depends on the size and lifetime on the collision fireball, as reflected by the clear dependence of 
the radial flow velocity on the impact parameter of the collision system in the right panel of Fig. 30. The 
larger fireballs formed in central collisions live longer, develop more radial flow and cool down to lower 
kinetic freeze-out temperatures than the smaller fireballs from more peripheral collisions. The chemical 
decoupling temperature T c h em , on the other hand, shows no dependence on the collision centrality (size 
of the fireball) at all; it can therefore not be controlled by a similar competition between scattering 
and expansion as seen in the kinetic decoupling temperature. T c h e m and Tki n rest conceptually on 
different footings. T c h em is not a kinetic decoupling temperature, but a Lagrange multiplier in the 
Maximum Entropy framework whose value indicates the critical energy density at which hadrons are 
formed. 

So if, as we concluded from the momentum spectra, after hadronization the system continues to 
expand and cool to lower temperatures, why doesn't the chemical decoupling temperature T c h em follow 
suit? The answer must be that, after being formed, the hadrons are unable to change their abundance 
ratios any further by inelastic rescattering. If this were not the correct interpretation, it would be very 
difficult to see how about a dozen different hadron species with abundances varying over 3 orders of 
magnitude could conspire to look like chemical equilibrium at T«T cr . This would be an unlikely 
accident indeed. In Sec. 6.4 we will see theoretical arguments why chemical decoupling has to happen 
more or less directly at T cr . 

6.3 Strangeness enhancement 

The fact that the chemical decoupling temperature has an uncanny resemblance with the predicted phase 
transition temperature is, in fact, only half the story. A similar picture of statistical hadronization at the 
critical energy density e cr arises even from an analysis of e + e~, pp and pp collisions [15]. What is really 
dramatically different in heavy-ion collisions is the level of strangeness saturation reflected in the appar- 
ent chemical equilibrium state: The left panel of Fig. 31 shows that the overall fraction of strange par- 
ticles is about twice as high in heavy-ion collisions as in elementary particle collisions ! In other words, 
strangeness is suppressed relative to its chemical equilibrium saturation value in elementary particle col- 
lisions, but this strangeness suppression has disappeared in ^4+^4 collisions. According to the preceding 
discussion this extra strangeness cannot have been produced by final state hadronic rescattering; it thus 
reflects the properties of the prehadronic state. This points to a new, fast strangeness production mech- 
anism, either before or during hadronization. Rapid strangeness production from thermal gluons was 
predicted as one of the key characteristics of a QGP [16]. 

A large fraction of the strangeness suppression seen in elementary particle collisions can actu- 
ally be understood within the statistical approach. Strong interactions conserve net strangeness exactly, 
i.e. S = without any fluctuations. This means that in principle strangeness conservation should not 
be treated in the grand canonical ensemble as in Eqs. (84) and (85), where it is conserved only on av- 
erage, but in the canonical ensemble were S is fixed exactly. In the infinite volume limit this makes 
no difference, but for the small "fireball" volumes created in elementary particle collisions it is a major 
correction. Strange hadron get suppressed in this way since, together with the strange hadron, always a 
second hadron with balancing strangeness has to be created inside the same small volume at the same 
time, and this requires more energy. However, even if this canonical effect were the complete explanation 
of the suppression in e + e~ and pp collisions (which seems quantitatively unlikely [152]), the disappear- 
ance of the effect in heavy-ion collisions would still be most interesting: The absence of canonical 
strangeness suppression in nuclear collisions tells us that creation of a strange hadron at a given position 
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Fig. 31: Left: The strangeness suppression factor of produced strange vs. nonstrange quarks, A s = 2 (ss) / (uu-\-dd) | produced, 
in elementary particle and heavy-ion collisions as a function of ^J~s [145]. The two points each for pp collisions reflect the 
inclusion (exclusion) of the initial valence quarks. Right: Enhancement factor for the midrapidity yields per participating 
nucleon in 158 A Ge V/c Pb+Pb relative to p+Pb collisions for various strange and non-strange hadron species [151]. 



in the fireball does not require the production of a particle with balancing strangeness nearby (as in pp), 
but strangeness can be balanced by production of an anti-strange hadron on the other side of the nuclear 
fireball! Now we know that all microscopic QCD processes which create ss pairs are local processes, i.e. 
the pairs are created at the same point. If the final strange hadron abundaces don't know about this any- 
more, but behave as if the strange and antistrange hadrons were created independently and statistically 
distributed over the entire nuclear fireball, the necessary implication is that they have lost their memory 
of the locality of the primary QCD process from which they came and have been able to communicate 
over macroscopic fireball distances. I cannot imagine how this is possible without a significant amount 
of strangeness diffusion before hadronization. 

A very striking way of plotting the "strangeness unsuppression" is shown in the right panel of 
Fig. 31 [151]: relative to p+Pb collisions, the number of produced strange hadrons per participating nu- 
cleon is the more strongly enhanced the more strange (anti)quarks for its formation are required. For 
and Q this enhancement factor is about 15! The tendency shown in Fig. 31 is completely counterintu- 
itive for hadronic rescattering mechanisms, where multistrange (anti)baryons are suppressed by higher 
thresholds than kaons and A's; but it is perfectly consistent with a statistical hadronization picture [153] 
where multi-strange particles profit more from the global strangeness enhancement than singly strange 
hadrons. 



6.4 Chemical kinetics after hadronization 

The kinetic freeze-out temperatures of about 120MeV extracted from he hadron momentum spectra tell 
us that after hadronization at T cr « 170 MeV the hadrons continue to rescatter for quite a while. This 
rescattering can be simulated in a hadronic cascade such as URQMD [154], using everything that is 
known about hadron masses and cross sections from the Particle Data Tables. These tables tell us that 
hadron cross sections are strongly dominated by resonances at well-defined collision energies, for ex- 
ample 7r+ N — » A — > tt+N, 7T+7T — ► p — ► 7r+7r, tt+K — ► K* — > 7T+K, etc. These resonances have large 
cross sections, hence these types of collisions happen frequently and are able to keep up with the col- 
lective expansion of the fireball by continually re-equilibrating the hadron momentum distributions to 
the falling temperature. However, almost all of these resonances have the tendency to decay again into 
exactly the same hadrons from which they were created (up to charge exchange). The cross sections 
for strangeness exchange processes are considerably smaller, and even smaller are inelastic hadronic 
processes which create or annihilate strange quark-antiquark pairs. The resonant collisions thus do not 
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Fig. 32: Time-dependence of the chemical compo- 
sition of an expanding hadronic fireball which was 
initialized at T CI with chemical equilibrium abun- 
dances [154], Unstable hadrons are seen to decay 
as stable hadrons approach their asymptotic abun- 
dances (top). In the bottom panel, the squares in- 
dicate the initial chemical equilibrium yields for 
stable hadrons, after taking into account feeddown 
from unstable resonance decays, whereas the cir- 
cles indicate the final stable hadron yields including 
hadronic rescattering. 



modify the total measured yields of pions, kaons, and nucleons. For example the total pion number 



N, 



total 



N w + 2N p + N A + N K * + ... 



(86) 



is frozen in at T cr , and similarly for other stable hadrons. This is seen in the bottom part of Fig. 32 which 
shows that, if URQMD is initialized with thermodynamic equilibrium abundances and a hydrodynmiccal 
generated transverse flow profile similar to the discussion in Sec. 4, the stable particle yields hardly 
change if they are calculated directly at T cr by simply letting all unstable resonances decay, but switching 
off all collisions, or if they are calculated at the end of the entire hadronic cascade. 

Looking a bit more carefully at the bottom part of Fig. 32 we see that URQMD tends to lose pre- 
dominantly (anti)baryons (and thus a fraction of the initial enhancement of multi-strange (anti)baryons) 
by baryon-antibaryon annihilation during the rescattering stage. It was recently shown [155] that this is 
to a large extent a manifestation of the lack of detailed balance in the codes which include processes like 
pp — > mr (with n = 5 — 6) but not their inverse. Rapp and Shuryak [155] argue that, as the system 
cools below T cn em> pions and kaons don't annihilate but instead build up a positive chemical potential 
which enhances the probability for the inverse reaction and strongly reduces the net annihilation of an- 
tibaryons. This is really fortunate, because it is this lack of abundance-changing processes during the 
hadronic expansion stage which allows us to glimpse the hadronization process itself through the final 
hadronic abundances, in spite of intense, resonance-mediated elastic rescattering among the hadrons 
between hadronization at T c h em ~ 170 MeV and kinetic freeze-out at Tf « 120MeV. 

Of course, the abundances of unstable resonances are not frozen in at T cr (see top panel of Fig. 32): 
Due to their strong coupling to the cooling pion fluid (to which cooling they actively contribute), their 
abundances readjust to the decreasing temperature. Note, however, that their abundances don't fall with 
the usual free-space exponential decay law, but more slowly since they keep being recreated by resonant 
hadronic interactions. But as the temperature drops, their abundances decrease. Detailed balance dic- 
tates that, the larger their pionic decay width, the later the resonance yields decouple and the smaller a 
temperature their decoupling abundances reflect. If we can reconstruct the resonances from their decay 



products we therefore expect resonance yields which are smaller than predicted from the thermal fit of 
the stable hadron yields: 

'K*\ fK*\ 
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Such studies are presently under way. 



6.5 Primordial hadrosynthesis: measuring the critical temperature for deconfinement 

We can summarize our insights gained in this chapter by stating as our Fourth Lesson: The hadron 
abundances "freeze out" directly at hadronization. What we see in the SPS and RHIC data are hadron 
abundances established during primordial hadrosynthesis via a statistical hadronization process. The 
measured abundance ratios give for the hadronization temperature Th a( j T cr 170 Me V, confirming 
predictions from lattice QCD. The QCD phase transition temperature has thus been measured. This 
measurement was possible because at SPS and RHIC the fireball expands so rapidly that in the late 
hadronic rescattering stage no inelastic chemical processes happen anymore. In this way the hadron 
abundances open a window onto the hadronization phase transition even though their momentum distri- 
butions continue to change and cool for several more fm/c. 



Fig. 33: QCD phase diagram with chemical 
freeze-out points extracted from hadron abun- 
dance ratios from heavy ion (A — 200) colli- 
sions at center of mass energies -^/snn rang- 
ing from 2.4 to 200 GeV [54]. The curve and 
shaded region indicate estimates for the decon- 
finement phase transition at finite baryon chem- 
ical potential /is from recent lattice QCD cal- 
culations [2, 156]. The critical endpoint where 
the transition changes from first order (high /is) 
to continuous crossover (low fis) is also indi- 
cated (although for unrealistically large quark 
masses [2]). Figure taken from Ref. [54]. 
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Figure 33 shows the chemical freeze-out points extracted from hadron yields in A+A collisions 
with A ~ 200, compiled for heavy-ion fixed-target and collider experiments from ^/snn ~ 2.4 GeV (SIS) 
to yj snn ~ 200 GeV (RHIC), and compares them to the latest estimates for the hadronization phase 
transition from lattice QCD [2, 156]. While for the SPS and RHIC the chemical decoupling temperatures 
are consistent with the predicted hadronization temperatures, the values extracted at the Brookhaven 
AGS and in particular at the SIS at GSI lie considerably below the phase transition. In these cases 
statistical hadronization can not explain the values for T c h em - H is likely that, due to the larger net baryon 
densities in lower-energy heavy-ion collisions, inelastic hadronic rescattering processes happen faster 
than at higher energies and are able to lead to kinetic readjustment of the chemical temperatue below 
T cr . One should also note, however, that the "quality" of the apparent chemical equilibrium seen at the 
AGS and SIS is not as good as at the higher energies since many of the heavier hadrons (in particular the 
strange and multistrange antibaryons) are too rare to be measured reliably. The data points labelled by 
"SIS, AGS" in Figure 33 thus comprise a much smaller number of measured hadron yields, in particular 
in the strange sector where only K, K, A and ip are available at these energies. 




7 PARTON ENERGY LOSS AND JET QUENCHING 



In Section 4 we saw that in Au+Au collisions at RHIC the soft momentum particles with p± < 2 GeV/c 
behave hydrodynamically. For larger pj_ the hydrodynamic description gradually breaks down. We saw 
this in Sec. 4.2 when discussing the elliptic flow and how it breaks away from the hydrodynamically pre- 
dicted increase with p± above p± ~ 2 GeV/c. In the single particle spectra the transition from collective 
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Fig. 34: Charged hadron transverse mo- 
mentum spectra in p+p and Au+Au col- 
lisions at ^/s^ = 200 GeV [157]. The 
Au+Au data are binned in different cen- 
trality classes as indicated. For central 
and semicentral Au+Au collisions, note 
the transition from an exponential be- 
haviour at p± < 2 GeV/c to a power law 
for p± > 4 GeV/c. 



hydrodynamic behaviour to hard scattering manifests itself by a change of shape: above p± > 3 — 4 GeV/c 
the spectra change from a thermal exponential shape to a power law as predicted by perturbative QCD. 
Figure 34 shows this for Au+Au collisions at RHIC, together with a comparison spectrum from p+p 
collisions. This power law affects fewer than 0.1% of all produced hadrons, but it is these rare high-p^ 
hadrons on which we will focus our attention in this section. 

These high-p^ hadrons stem from the fragmentation of even-higher-p^ partons which, as noted 
in Sec. 3.1, are created very early in the collision, at Tf orm — l/p±- On their way out of the evolving 
collision fireball they interact with the hot matter and probe its properties (density, opacity, etc.). In this 
respect they are very similar to positrons in PET (= Positron Emission Tomography) where one injects a 
positron-emitting source into some living organ and then explores properties of that organ by detecting 
these positrons and their energies outside the body. Measuring the properties of hard hadrons emitted 
from a relativistic heavy-ion collisions thus corresponds to Parton Emission Tomography, but since this 
would have the same acronym which might lead to all kinds of confusion, and since we don't really 
detect partons but rather their hadronic fragments which form a jet, I follow M. Gyulassy [158] and call 
this method JET (for Jet Emission Tomography). 

7.1 Radiative energy loss of a fast parton 

The idea of jet quenching by parton energy loss goes back to an unpublished preprint by J.D. Bjorken in 
1982 [159], but he incorrectly identified elastic parton scattering as the dominant medium interaction by 
which the fast parton loses energy. The correct mechanism, namely induced gluon bremsstrahlung, was 
identified by Gyulassy and collaborators [160] and first fully evaluated by Baier, Dokshitzer, Mueller, 
Peigne, and Schiff [161]. Its quantitative effect on the energy loss has now been calculated by several 
groups under various approximations (see [162] for a recent review). 

In QCD the dominant process by which a fast parton with E paiton S> 1 GeV loses energy is by in- 
duced gluon radiation, see Fig. 35. The color charge of the fast parton (which could be a quark or a gluon) 
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Fig. 35: Typical gluon radiation diagram (from [162]). 



interacts with the color charges of the medium (which can be modeled as static external charges since 
their motion can be neglected relative to that of the fast particle) and, as a result, emits a bremsstrahlung 
gluon. What makes things different from and more complicated than the analogous process in QED is 
the fact that the emitted gluon itself carries color charge and again interacts with the color charges in the 
medium (see Fig. 35), and that all these medium interactions must be added up coherently. In contrast 
to the fast particle, whose trajectory can be taken as a straight line, the reinteractions of the much softer 
emitted gluon with the medium induce a random walk in its transverse momentum fcx- The reinterac- 
tions of the emitted gluon with the medium are characterized by a mean free path A g , also indicated in 
the Figure. 

In the limit E pSLrton — ► oo and for a thick target, L » A g , BDMS [161] found for the energy loss 



where Cr is the color Casimir eigenvalue of the fast parton (= N c for gluons), fi is the color Debye 
screening length of the medium (controlled by its density), and v ~ 1—3 is a factor which depends loga- 
rithmically on the thickness L of the medium and the gluon mean free path, v ~ ln(L/A g ). The ratio 



plays the role of a momentum transport coefficient of the medium and exhibits the explicit dependence 
of the energy loss on the density p of the medium. The non-Abelian nature of QCD reflects itself in the 
characteristic non-linear dependence of the energy loss (88) on the thickness L of the medium. The extra 
factor of L compared to the naive expectation arises from the random walk in of the emitted gluon 
which increases its probability to separate itself (decohere) from the fast parton linearly with L. 

The asymptotic expression (88) was improved for "thin" plasmas and finite energy kinematic 
effects by Gyulassy, Levai and Vitev (GLV) [163] by use of an expansion in powers of the "opacity" of 
the medium (i.e. the integral over density times cross section along the path of the fast parton, n = Lj A g ). 
Their leading order contribution, which strongly dominates the energy loss, is 

While the characteristic quadratic L dependence of Eq. (88) is not altered by keeping only the first term 
in the opacity expansion, the dependence of the energy loss on the parton energy E = E'parton is quite 
different, due to the factor N(E) in the denominator. For E — > oo this factor approaches the BDMS 
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(88) 




(89) 



value, N(E) — ► 4, but as the parton energy E decreases N{E) increases rapidly, suppressing energy loss 
for low-E partons. In the range 2 GeV < E < 10 GeV, for fixed L the fractional energy loss AE^ V /E 
is almost constant [164]. 



7.2 Soft vs. hard particle production 

According to QCD and asymptotic freedom, the production of high-p^ partons can be calculated per- 
turbatively and is proportional to the number of binary nucleon-nucleon collisions which for a nucleus- 
nucleus collision A+B at impact parameter b is given by 

JVcoii(&) = dxdyT A (x+b/2,y) ■ T B (x-b/2,y). (91) 

Here ao is the total inelastic nucleon-nucleon cross section, the integral is over the transverse plane, and 
Ta(x, y)= J Pa(x, y, z) dz is the nuclear thickness function of a nucleus of mass A with density pro- 
file pa- The impact parameter b can be estimated from the total charge multiplicity, by equating equal 
fractions of the total multiplicity divided by the maximum multiplicity measured in the most central col- 
lisions with the corresponding fractions of the total cross section, calculated geometrically as a function 
of b [165]. Due to the fact that for fixed impact parameter B the number of produced hadrons fluctuates, 
this determination of b from dN^/dy is uncertain by a fraction of 1 fm, which induces a corresponding 
uncertainty in the relation between iV co n and dN^/dy [165]. 

The argument that hard particle production should scale ~ N co \\ exploits the fact that hard particles 
are produced on short time scales r ~ l/p±, and that hard particle production on successive nucleons in 
the nucleus therefore happens incoherently. The same is not true for soft hadron production which 
involves the coherent scattering of a projectile nucleon with several target nucleons, due to the Landau- 
Pomeranchuk-Migdal (LPM) effect of a finite formation (or "decoherence") time (see the last paper in 
[160] for an explanation of the LPM effect). The net result of the LPM effect is a reduction of soft particle 
production by destructive interference, resulting in its phenomenologically observed approximate scaling 
with the number N paxt of participating (or "wounded") nucleons: each struck nucleon contributes only 
once to soft particle production, and suffering more than one collision in sequence does not increase the 
soft particle yield. The number Ap art of wounded nucleons can again be calculated geometrically, by 
a formula similar to Eq. (91) but involving the sum of the nuclear thickness functions instead of their 
product (see e.g. Ref. [63]). 



7.3 Suppression of high-pj^ hadrons in central Au+Au collisions at RHIC 

In the 5% most central Au+Au collisions at v / snn = 200 GeV at RHIC, the average number of partici- 
pating nucleons is (N part ) = 344 (i.e. 172 times higher than in a p+p collision) whereas the number of 
binary nucleon-nucleon collisions is (AT coll ) = 1074 (i.e.1074 times higher than in a p+p collision) [165]. 
The difference between "participant scaling" (expected for soft particles) and "binary collision scaling" 
(expected for hard processes) is in this case a factor of 1074/172 = 6.2. The "nuclear modification factor" 

J_ dN AA /dydp ± 
~ Kifi dNPP/dydp ± 

for A = 197 (gold) should therefore be Raa = 1 for particles whose production scales with N co \\, and 
Raa — 0.16 for particles whose production scales with Apart. Hence, as a function of increasing p±, we 
should expect Raa to rise from a value near 1/6 at low p± to about 1 at high p±. 

However, it has been known for over 25 years from minimum bias p+^4 collisions at Fermilab 
[166] that Raa reaches values even larger than 1 at p± > 2GeV/c. The explanation of this "Cronin 
effect" goes as follows: At high transverse momenta the p± -spectrum of hadrons produced in p+A is 
not only shifted up in normalization by a factor iVcoii ~ A as appropriate for a hard scattering process, 



but on top of this also horizontally towards higher p± (which at fixed p±, of course, manifests itself as 
an additional gain, due to the falling nature of the p^-spectrum). The horizontal shift towards larger p± 
arises from the fact that the partons inside the projectile nucleon, before making the hard collision which 
ends up producing the measured high-p^ hadron, already suffer multiple elastic collisions with other 
target nucleons which they encounter first, thereby acquiring transverse momentum k± which grows in 
a random walk with the square root of the number of these elastic collisions. Once the hard inelastic 
collision happens, the projectile parton already brings in this "initial k±", thereby giving and extra k± 
kick to the produced hard parton. Since this extra k± is a decreasing fraction of the observed p± as 
p± becomes larger, the Cronin enhancement should disappear as p± — > oo. (For the same reason the 
Cronin effect should become weaker as ^/snn increases.) Hence, for large p±, Raa should indeed 
approach the value 1, but not from below due to soft scaling at low p±, but rather from above due to 
Cronin enhancement at intermediate p± [167]. This behavior is seen in the two top curves (labelled "no 
dE/dx") in the right part of Fig. 36. 
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Fig. 36: Left: The nuclear modification ratio Raa defined in Eq. (92) for neutral pions measured in a+a collisions at the 
ISR [168], in Pb+Pb collisions at the SPS [169], and in Au+Au collisions at RHIC [170, 171]. (Figure taken from Ref. [172].) 
Right: Raa for ^° in central Au+Au collisions at RHIC [170] together with theoretical calculations [164, 167, 173]. (Figure 
taken from Ref. [174].) 



The left panel of Figure 36 shows the measured nuclear modification factor for neutral pion pro- 
duction in Pb+Pb collisions at the SPS and in Au+Au collisions at RHIC. For comparison also results 
from a+a collisions at the ISR in the range 4GeV/c<p^ < 7GeV/c are included which exhibit the 
Cronin effect in this p± region. The Pb+Pb data at y / SNN = 17 GeV show a dramatic growth of Raa 
with increasing p±, starting from participant scaling at low p± but then exceeding the binary collision 
limit and exhibiting a strong Cronin enhancement above pj_ > 2 GeV/c. (This Cronin effect is larger than 
in the ISR data, due to both the larger collision system, implying a larger initial k±, and the lower center 
of mass eneergy.) 

In stark contrast, the RHIC Au+Au data at both 130 and 200 GeV cm. energy hardly rise at all 
above the participant scaling level. With some good will R\ A is seen to feature a small peak around 
p± simeq 2 GeV Raa, but at a level of less than 50% of the binary collision limit. At higher p±, instead 
of increasing towards 1, Raa decreases again towards a value of about 0.2. No sign of the Cronin ef- 
fect and of binary collision scaling anywhere! If we jump ahead to Fig. 37 which shows in the lower 
parts of both panels (labelled by "Au+Au") the p± -dependence of Raa for all charged particles, and 
if we assume that charged and neutral pions (triangles in the left panel) have the same nuclear modifi- 
cation factor, we see that Raa for heavier charged hadrons must be larger than for pions in the region 
1 GeV/c < p± < 4 GeV/c, but also drop back to well below 1 for p± > 4— 5 GeV/c. 

The conclusion we must draw from these data is that in central Au+Au collisions at RHIC high-p^ 
hadron production is suppressed by at least a factor 5 compared to the perturbative QCD prediction of 



binary collision scaling plus Cronin effect (top curves in the right panel of Fig. 36). In fact, looking at 
various collision centralities in Au+Au collisions, it was shown explicitly by the PHOBOS Collaboration 
[175] that above p± > 4GeV the charged hadron spectra scale again linearly with A^ part (i.e. like a soft 
process), instead of scaling with N co \\ as expected for a hard process. Furthermore, a crucial control 
experiment has recently been completed where the nuclear modification factor was measured in d+Au 
collisions (top parts of Fig. 37). One sees that in deuteron-gold collisions Raa = -^dAu increases from 
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Fig. 37: The nuclear modification ratio Raa for neutral pions (left) and charged hadrons (left and right), measured in d+Au 
and central Au+Au collisions at v /55^= 200 GeV, by PHENIX [176] (left) and by STAR [177] (right). At intermediate pj_, 
1.5GeV/c<p_L <4GeV/c, the Raa is systematically larger for charged particles than for it , indicating a larger Raa for 
kaons and protons than for pions in this px-range. 

participant scaling at low p± to above 1 at p± ~ 1.5 GeV/c and shows the expected Cronin enhancement 
(rather than a suppression) which persists at least up to p± = 8 GeV/c. Again all charged particles taken 
together (which includes heavier hadrons such as kaons and protons) show a larger Cronin enhancement 
than pions alone (left panel in Fig. 37), and in central d+Au collisions the Cronin effect is a bit larger 
than in minimum bias collisions (right panel in Fig. 37), as expected from the random walk argument for 
the initial k±. 

So the suppression of high-p^ hadrons in central Au+Au collisions must be a final state effect 
which involves the dense matter created in the heavy-ion collision. It can not be explained as an initial 
state effect, i.e. by some property of the internal structure of the gold nucleus which shows up when 
probing it at high energies (such as the existence of a "color glass condensate", see the short discussion 
in Sec. 3.4 and the recent review [28]): such an effect should still be visible (i.e. some suppression should 
persist) if only one of the colliding nuclei is gold [178]. 

The right panel in Figure 36 shows that only theories which include an energy loss dE/dx for 
the fast partons that produce the high-p^ hadrons can account for the Raa data at RHIC. Although the 
different theoretical predictions shown in the graph, which use different approximations for the imple- 
mentation of the radiative energy loss, give somewhat different p± -dependences of Raa and do not agree 
perfectly with the data, it is obvious that a significant amount of energy loss is required to account for 
the observed suppression by more than a factor 5 relative to the expected binary collision scaling plus 
Cronin (upper curves). 

The observed scaling with A^ part instead of iV co ii suggests a very interesting geometric interpre- 
tation: Due to the approximately constant density inside an atomic nucleus, the number of participating 
nucleons is proportional to the volume of nuclear matter affected by the collision, Np art oc Vg^ all . The 
number of binary collisions scales very accurately as iV co ii oc N^ t . So the difference between partici- 
pant and binary collision scaling is a factor N pi ^J t 3 oc oc l/i?fi re b a ii> ie. a surface/volume ratio: 
the high-px hadron production data in central Au+Au collisions at RHIC are consistent with the sim- 



pie hypothesis that high-p^ hadrons are only emitted from the surface, but not from the interior of the 
fireball! Now why should this be the case? 

7.4 Jet quenching in central Au+Au collisions at RHIC — "JET of the QGP" 

Experiments provide a unique answer also to this question [179]: Jets from partons which have to travel 
a significant distance through the dense fireball matter formed in the collision lose so much energy that 
they are no longer recognizable as jets and become part of the hot matter of soft particles. To show 
this experimentally one must first find jets. This is not easy due to the huge number of soft hadrons with 
p± < 2 GeV which contribute so much transverse energy that the usual method of finding a spike in Ei in 
a certain small bin of solid angle fails. Another method for finding jets does work, however: One triggers 
on a fast particle (say, with 4 < p± < 6 GeV/c [179]) and then looks for azimuthal angular correlations 
with other, not too soft hadrons (say, with p± > 2 GeV/c, in order to remove most of the background 
from uncorrected soft hadrons). Since fast partons fragment into a jet of hadrons pointing within a 
relatively narrow angular cone ("jet cone"), jets manifest themselves by positive angular correlations at 
small angles relative to the fast trigger particle. This is shown in Fig. 38 where it gives rise to the peak at 
A(p = 0. 
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Fig. 38: Jet suppression in central Au+Au collisions at ^snn = 200 GeV at RHIC. Left: Azimuthal angular correlations in 
Au+Au collisions, in four centrality bins, compared with properly scaled p+p data to which the effect of elliptic flow V2 
measured in the same p±- and centrality window has been added. (For details see Ref. [179].) Whereas such a superposition 
of p+p data and elliptic flow can reproduce the peripheral (60-80%) Au+Au data, it overpredicts the far-side correlations at 
Atp — n for semicentral and central collisions. Right: Azimuthal angular correlations in p+p as well as min. bias and central 
d+Au collisions (upper diagram), and in p+p, central d+Au and central Au+Au collisions (lower diagram) [177]. For details 
about the subtraction procedure in the bottom part see the original paper [177]. 

In perturbative QCD, hard partons are produced in pairs with 180° opening angle in the pair 
center of mass frame. This leads to a second azimuthal angular correlation peak at Acp = n. The solid 
histogram in the right panel of Fig. 38 shows this for p+p collisions at ^/snn = 200 GeV [177]. The "far- 
side" peak is flatter and wider than the "near-side" peak, due to well-understood trigger effects and a 
non-zero probability for semihard gluon emission by one of the two fast partons. Except for a somewhat 
higher background from uncorrected pairs due to increased soft particle production ("pedestal effect"), 
the same pattern is seen in minimum bias and central d+Au collisions (Fig. 38, upper part of the right 
panel). 



However, as one goes to Au+Au collisions and moves from peripheral to central collisions (bottom 
to top in the left panel of Fig. 38), the away-side peak at Atp = it disappears. Whatever enhanced angu- 
lar correlation survives in almost central Au+Au collisions at Atp = it is perfectly compatible with the 
measured elliptic collective flow V2 in the p± -range where the angular correlations are measured [179]. 
Whereas the far-side jet correlations should be localized in a narrow window of relative pseudorapidity 
Arj between the correlated particles, elliptic flow is a global collective effect which is correlated with 
the reaction plane and therefore present at both small and large relative pseudorapidity. By subtracting 
the correlations for pairs with Arj > 0.5 from those for pairs with Arj < 0.5 one in fact completely re- 
moves the elliptic flow contribution [179]. For semicentral Au+Au collisions the result of this procedure 
is a strongly reduced far-side angular correlation peak, and in central Au+Au collisions it is completely 
absent (lower diagram in the right panel of Fig. 38)! 
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Fig. 39: Left: Suppression of near-side (squares) and far-side angular correlations (circles) in Au+Au collisions at RHIC, 
for two different lower p±_ cuts for the trigger particle [179]. Note the disappearance of the far-side correlations in central 
collisions (i.e. for large numbers iVp ar t of participating nucleons). Right: Cartoon illustrating outward jet emission from the 
fireball surface and quenching of its inward-moving partner by the dense interior. 

One can further quantify this by comparing the area under both near- and far-side peaks of the 
azimuthal correlation function (after subtracting elliptic flow) with the corresponding areas in p+p colli- 
sions. Figure 39 shows the ratio of these areas for the near-side peak as squares and for the far-side peak 
as circles. For very peripheral Au+Au collisions both ratios are close to, but slightly less than 1. The 
reduction below 1 can be understood as a trigger effect [180]. As the Au+Au collisions become more 
central (i.e. -/V par t increases), the near-side correlations seem to increase, in particular for lower values 
of the trigger-px, which may indicate an incomplete subtraction of elliptic flow effects. The away-side 
correlations, on the other hand, get weaker in less peripheral collisions and completely vanish in central 
collisions. If the elliptic flow were larger than assumed in the subtraction, this would bring down both 
the near- and far-side correlation strength, rendering the disappearance of the far-side jets even more 
dramatic. 

The interpretation suggested by these data is indicated by the cartoon in the right part of Figure 39: 
The high-px trigger hadron identifies a jet moving outward from the fireball. The yield of these high-p^ 
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hadrons scales by a factor N p J t oc 1/ -Rgrebaii more weakly than expected from perturbation theory, so 
these outward-moving jets come only from a thin surface layer. Fast partons created in the interior and 
trying to move out, as well as the inward-moving partners of the observed outward-moving jets, lose 
so much energy in the dense medium that they no longer make hadrons with p±>2 GeV which are 
included in the angular correlation function. (If they still produced hadrons with p± > 2 GeV and had 
lost no energy but only their angular correlation with the fast trigger particle, the total yield at high p± 
should still scale with N co n, i.e. Raa should be 1.) 



This amounts to "complete jet suppression": The hot and dense medium (which we know from 
Sec. 4 must be a thermalized QGP) suppresses all jets from fast partons except for a small fraction 
being produced near the surface and moving outward. Since in jet fragmentation the leading hadron 
carries on average about half the p± of the hard parton, a trigger hadron with a pj_ between 4 and 6 
GeV/c corresponds to p± ~ 10 for both primary partons. The complete disappearance of the far-side jet 
in central Au+Au collisions requires the far-side parton to lose at least 6 GeV of momentum such that 
its fragments are all below the 2 GeV/c cut imposed on the angular correlation function. The energy 
lost by these partons thus becomes part of the thermalized and hydrodynamically expanding medium at 
p± <2GeV/c. 

So here is our Fifth Lesson: The quark-gluon plasma strongly suppresses jets. Inward-moving 
partons with p± < 10 GeV don't make it to the opposite edge and become part of the low-p± "thermal 
soup". Only outward-moving jets formed near the fireball surface survive. The era of JET (Jet Emission 
Tomography) of the quark-gluon plasma has begun: by investigating the A-, b- and emission-angle 
dependence of the far-side jet suppression and combining it with information extracted from soft particle 
emission we should be able to accurately determine the fireball conditions (density and temperature) in 
the early stage. 

8 CONCLUSIONS AND DISCLAIMERS 

The title of these lectures has been "Concepts in Heavy-Ion Physics", but these lecture notes do not 
give a comprehensive overview of all the relevant concepts. In 3 x 75 minutes one simply can't cover 
everything. Notable omissions include (I give references to good recent reviews) J/ijj suppression [181] 
as well as direct photon and thermal dilepton emission [13, 182]. While RHIC has not yet produced 
sufficient data on the first two of these items, J/ip suppression has been a hot subject at the SPS, and 
some very interesting signals were also seen in the photon and dilepton spectra [10, 183]. I also omitted a 
more detailed discussion of the very interesting quark-coalescence approach to the production of hadrons 
at intermediate transverse momenta (1.5 GeV /c < p± < 6GeV/c) where some really nice systematics has 
been discovered which strongly indicates that in this kinematic region quark counting rules, causing a 
meson-baryon splitting of the elliptic flow and hadron yields instead of the hydrodynamic mass splitting 
observed at lower p± [184]. I may include a more detailed discussion of these topics in some future 
lecture notes. 

Relativistic heavy-ion physics is making progress in huge strides. These lectures have taught us 
a few important lessons which, I think, are now solidly established, but with the steady stream of RHIC 
data new, smaller lessons are being taught to us by Nature every week and keep complementing the 
picture literally as I sit here writing these notes. (Just compare these notes and the many new data I 
present here with the copied transparencies distributed at the summer school! And this write-up doesn't 
even mention everything that happened during the last year . . . ) Nevertheless, let me repeat once more 
what I think are the important and firm conclusions one can draw at this point in time (note that, in spite 
of the new evidence and significantly improved data, these still agree with what I said at the school over 
a year ago): 

• The heavy-ion reaction zone undergoes violent explosion; this proves the existence of strong 
thermal pressure and intense rescattering of the quanta produced in the collision. The bulk of 
the hadron data at p± < 2 GeV is well described by hydrodynamics, and decoupling occurs at 
Tf rj 100 - 120 MeV with average transverse flow (/3 T ) rj 0.5-0.6 (=> "The Little Bang"). 

• The measured elliptic flow in non-central collisions almost exhausts the hydrodynamic upper limit. 
This provides stringent limits on the time scale for thermalization, r t h C rm < 1 fm/c, and shows that 
the quark-gluon plasma is a strongly coupled fluid and thermalization in the QGP is controlled by 
non-perturbative processes. 

• At T t hcrm the energy density of the thermalized system exceeds the critical value for quark decon- 
finement by more than an order of magnitude. Hence, this thermalized state must be a Quark-Gluon 



Plasma. At RHIC energies, the QGP lives for about 5— 7fm/c before hadronizing. At the time of 
hadron formation the elliptic flow has already saturated. 

• Soft hadronization at T cr i t ~ 170 MeV is a statistical process which generates a chemical equi- 
librium distribution. This equilibrium is not due to kinetic processes involving inelastic scat- 
tering between hadrons, but the result of a statistical formation process following the principle 
of maximum entropy. Due to rapid expansion and dearth of inelastic processes among hadrons, 
the hadronic chemical composition decouples immediately. Therefore the critical temperature for 
hadronization T cr = T c h C m can be extracted from the final hadron yields. One finds T cr m 170 MeV, 
confirming lattice QCD predictions. 

• The QGP induces strong energy loss for high-p^ partons. Partons with p± < 10 GeV which 
travel inward through the fireball lose enough energy to become indistinguishable from the soft 
hadron background ("complete" jet quenching). Only outward-moving surface jets are emitted, 
explaining why high-p^ hadron production (p± < 8 GeV) is found to scale oc iVp ar t instead of 
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oc N co \\ oc -/Vpart as perturbatively expected. 
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